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ABSTRACT 


Pilot-vehicle-display systems theory is applied to the 
analysis of proposed vertical situation displays for manual 
control in approach-to-landing of a STOL aircraft. The effects 
of display variables on pilot workload and on total closed-loop 
system performance was calculated using an optimal-control 
model for the human operator. 

The steep (7.5 deg) approach of an augmentor wing jet STOL 
aircraft was analyzed. Both random turbulence and mean-wind 
shears were considered. Linearized perturbation equations were 
used to describe longitudinal and lateral dynamics of the air- 
craft. The basic display configuration was one that abstracted 
the essential status information (including glide-slope and 
localizer errors) of an EADI display. Proposed flight director 
displays for both longitudinal and lateral control were also 
investigated. 

It was found that with the basic EADI-status displays 
Category II window specifications would be exceeded more than 
5% of the time unless unacceptable (or unachievable) levels of 
pilot workload were demanded. The proposed interim flight 
directors improved both performance and workload; they reduced 
the failure probability by about a factor of two or allowed a 
fixed success probability to be achieved with much less atten- 
tional demand. However, further improvements in performance 
via display design appeared unlikely, though greater reductions 
in workload may be possible. 
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A preliminary approach to designina display command infor- 
mation using the optimal-control model was also investigated. The 
procedure yielded a director law that resulted in substantially 
reduced workload. 
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INTRODUCTION 

One of the design goals for an aircraft instrument panel is 
to minimize display workload. That is, acceptable system per- 
formance should be achieved with a minimum of the pilot's atten- 
tion reouired by the flight displays so that the pilot may attend 
safely to the problems of vehicle manacrement other than continuous 
control. Sophisticated displays such as flight directors and 
pictorial displays are sometimes included in the instrument panel 
to minimize workload, especially in approach and landing. Recent 
developments in pilot-vehicle-displav systems theory [1-6] indi- 
cate the feasibility of determining analytically the relative 
workload associated with various displays for approach and 
landina. In particular, it is possible to assess the reduction 
in pilot workload that might be provided by flight-directors and 
other forms of display augmentation. These same techniaues can 
also be used to predict closed-loop measures for evaluatina over- 
all system performance with the various displays. 

The purpose of this program was to apply the pi lot-vehicle- 
display systems theory developed in [1-6] to the analysis of 
proposed vertical situation displays for manual control and 
monitoring in approach- to-landing of a STOL aircraft. The 
effects of display variables on pilot workload and on total 
closed-loop system performance was to be calculated using the 
aforementioned techniques. 

The steep (7.5 deg) approach of an Augmentor Wing Jet STOL Re- 
search Aircraft (AWJSRA) was analyzed. Both random turbulence and 

wind shears were considered. Linearized perturbation eouaticns 
were used to describe longitudinal and lateral dynamics of the 
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aircraft. The basic display configuration was one that abstracted 
the essential status information (including glide-slope and loca- 
lizer errors) of the STOLAND-EADI display [7] (hereafter, referred 
to as EADI-status or, simply, status display) . Proposed flight 
director displays for both longitudinal and lateral control were 
also investigated. Allocation of pilot attention, performance 
at the approach-window and pilot-workload were all investigated. 

In the course of this effort some modifications and exten- 
sions of the pilot-vehicle-display system model of [1-6] and of 
the corresponding computer programs were implemented. Multi-control 
problems were considered for the first time; the wind-shears that 
were analyzed were more general than those previously investigated 
[6] ; and the programs were made more amenable to the investigation 
of time-varying situations and to workload analysis. 

A concise summary of the main results of the analysis, with 
references to supporting data, is presented in the concluding 
chapter of this report. It is strongly recommended that the 
reader with insufficient time to go through the details of the 
analysis initially should refer to the conclusion directly. 

For convenience, the main conclusions are highlighted below. 

1. The EADI-status displays are inadequate for the 
approach task in that the Category II window 
specifications are exceeded more than 5% of the 
time unless unacceptable (or unachievable) levels 
of pilot workload are demanded. 
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2. The proposed interim-flight directors improve both 
performance and workload. On the average, the flight- 
directors reduce the probability of a missed approach 
by a factor of two or, as a corollary, for a fixed 
probability of success the workload with the direc- 
tors is less than half that without them. 

3. Further improvements in performance via display design 
are likely to be small (using the measurements currently 
available) , but greater reductions in workload may be 
possible . 

4. Non-zero mean winds, with shear-variation, degrade 
performance by producing mean errors and by increasing 
response variability. 

The organization of the report is as follows. The analysis 
procedures are reviewed in Section 2. The basic EADI-status 
configuration is analyzed in Section 3. In Section 4, proposed 
flight director laws are investigated . Section 5 contains a 
discussion and analysis of an approach to design of flight 
director laws that is based to the procedures used in this 
study. Only the longitudinal director is considered there. 

In Sections 3-5, the analysis is restricted to "steady-state" 
gust regulation at the decision-height; the response to wind 
shears is analyzed in Section 6. Concluding remarks are pre- 
sented in Section 7 and pertinent mathematical details are 
included in three Appendices. 
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ANALYSIS METHODS 

The procedures used in this study to analyze vertical 
situation displays for STOL approach have been developed, for 
the most part, over the past several years. The techniques 
and the basis for them (the optimal-control model for pilot- 
vehicle-display systems analysis) are described in detail in 
References 1-6. Here, we summarize the important features of 
the system model. Procedures for determining closed-loop sys- 
tem performance and indicators of pilot-workload are also 
outlined. 

The pilot-vehicle-display system is illustrated in 
Figure 1. The model includes representations of vehicle dyna- 
mics, environmental disturbances, the display system and a 
model for the pilot. The model for the pilot is the so-called 
"optimal-control model" [1-6] . 

System Model 


Dynamics 

System dynamics are approximated by the following linear 
state equation: 

x(t) = A x(t) + B u(t) + E • w(t) (1) 

with initial condition 


E{x(t )} = x ; E { (x ( t ) - X )(x(t ) - x )'} = X (2) 

— O — O —\j — o — o — o 


Preceding page blank 
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DISTURBANCES ^(t) 



MOTOR OBSERVATION 

NOISE NOISE 


HUMAN OPERATOR MODEL 


FIGURE 1. Structure of Pilot-Vehicle-Display System Model 
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where x(t) is the vector of system states, u(t) is the vector 
of the pilot's control inputs, and w(t) is a vector of linearly 
independent white noise disturbances. The system dynamics 
include vehicle dynamics f any "dynamic" filtering associated 
with the generation of input disturbances (see below) and, if 
necessary, sensor, actuator or display dynamics. In this 
study, the matrices A, B and E are constant, and correspond 
to a steep-approach flight condition. The constancy of A, B 
and E is not a necessary requirement of the procedure. Piece- 
wise-constant matrices are allowed and may be used to explore 
some time-varying situations; however, this increases computational 
costs significantly and makes the interpretation of results more 
difficult. 

Environmental Disturbances 

If external forcing functions are rational noise spectra 
of first order or higher, as is the case for most turbulence 
models, they are represented by white noise (w) passed through 
a linear filter. Then, "input states" are augmented to the 
system state and the filter dynamics form part of the A matrix. 

Disturbances such as constant winds or wind-shears are 
modeled, essentially, by adding non-zero mean components to w; 
any dynamics associated with these disturbances are accounted 
for in the manner described above. For this situation, the 
model for the human operator must be modified so as to account 
for time-varying adaptation to such disturbances. These modi- 
fications are discussed in Appendix A. The equations needed 
for the computation of the response of the overall system to 
a particular "sample" mean-disturbance are developed in 
Appendix B. 


7 



Report No. 2484 


Bolt Beranek and Newman Inc. 


Display Variables 

The display-variables are assumed to be linear combinations 
of the state and control variables and are given by the "display 
vector" : 

^(t) = C x(t) + D u(t) (3) 

The matrices C and D may be time-varying (piece-wise constant) 
to account for changes in the quantities being displayed or 
"observed". For example, in the STOLAND-EADI display [7] anaular 
glide path deviations are presented in the early part of the 
approach, whereas linear glidepath deviations are presented 
as the runway threshold is neared. 

As noted above, display dynamics are included in the A 
matrix. For the present study the only "display dynamics" 
considered were those associated with generating the "flight 
director" laws investigated in Chapters 4-6. 

Task Requirements 

Task requirements are stated in terms of "cost weightings" 
associated with various system variables in a quadratic cost 
functional of the form 

J (u) = E {£* Q £ + u ' R u + U' G u) (4) 

It is assumed that the pilot selects his control response to 
minimize the appropriate J. 
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The selection of cost functional weightings may, in general, 
be based on objective or subjective factors. For relatively 
simple, single-variable control situations, good approximations 
to experimental measurements have been obtained with a cost 
functional consisting simply of a weighted sum of system error 
variance plus control-rate variance [2]. The cost on control - 
rate represents, in part, a subjective penalty imposed by the 
controller on making rapid control motions. In addition, this 
term may account indirectly for physiological limitations on the 
pilot's bandwidth. 

For complex multi-input, multi-output tasks, the cost 
weightings can not be chosen in so simple a fashion. There are 
several ways one can proceed. Values for the weightings can be 
selected to keep mean-squared output levels within prescribed 
tolerances; they can be assigned via pilot questionnaire; they 
can be chosen by trade-off analysis; finally, they may be assigned 
such that the resulting optimal closed-loop system has certain 
"desirable" properties. 

For the analysis performed in this study, weighting coef- 
ficients were selected on the basis of maximum allowable deviations 
(or limits) for the various problem variables [8], A unit amount 
of cost was associated with a given variable when the magnitude 
of the "error" (i.e., deviation from trim) was equal to the nomi- 
nal limit. Thus, the weighting coefficient for each variable was 
computed simply as the inverse of the square of the corresponding 
limit. Values for the limits were determined partly from Category 
II window specifications, partly from physical limitations and 
partly from a knowledge of human preference and capabilities. 

This method for choosing cost-weightings represents a departure 
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from our previous approach {e.g., see [6]). We believe that it 
has the advantage of being relatively straightforward to apply. 

A possible disadvantage is that it may involve selection of more 
weightings than are actually necessary to give a reliable predic- 
tion of performance. Finally, it is worth noting that it is 
possible to allow the weightings to vary with range to represent 
the pilot's tendency or desire to "tighten" his control as 
touchdown is approached (see Chapter 6). 

Pilot Model 

For purposes of this discussion it is convenient to 
consider the model for the pilot as being comprised of three 
parts : 

(i) An "equivalent" perceptual model that translates 
displayed variables ^ into delayed, "noisy" per- 
ceived variables y via the relation 

‘‘-P 

£ p (t) = £(t-T) + v y (t-T) (5) 

where x is an "equivalent" perceptual delay and 
v y is an "equivalent" observation noise vector.* 

(ii) An estimation and control-command generation pro- 
cess that consists essentially of a Kalman filter, 
a least mean-squared predictor and a set of 
"optimal gains". This models human information 
processing and compensation behavior. 

*The use of the word equivalent in this context is to emphasize 
that the parameters may be lumped representations of a variety of 
limitations that can not be "identified" separately by existing 
measurement techniques. 
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(iii) An equivalent "motor" or output model that accounts 
for possible "bandwidth" limitations of the human 
and his inability to generate perfect control res- 
ponses, by transforming "commanded" controls, u , 

— c 

into control inputs via the transformation 


T„ u + u 
— N — — 


U + V 

— c -m 


( 6 ) 


where T„ is an "equivalent" "Neuro-motor" lag-matrix 

and v is an "ecrui valent" motor-noise vector. 

— m 


The optimal predictor, optimal estimator, and optimal 
gain matrix represent the set of "adjustments" or "adaptations" 
by which the pilot tries to optimize his behavior. The general 
expressions for these model elements are determined by well- 
defined mathematical rules that are described in Reference 2. 


The time-delay, observation-noise, motor-noise, and, in 
some cases, the elements of the T N matrix represent inherent human 
limitations, and those of the human/system interface, that tend 
to constrain the range of the pilot's behavior. These limita- 
tions are discussed below in more detail. 


(a) Time Delay. The various internal time delays associ- 
ated with visual, central processing and neuro-motor pathways are 
combined and conveniently represented by a lumped equivalent per- 
ceptual time delay x. Typical values for this delay are 0.2 + 

. 05 sec. t 5] . 
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(b) Neuro-Motor Dynamics. We do not include "neuro-motor" 
dynamics directly among the inherent limitations of the human. The 
matrix of Equation (6) limits the "bandwith" of the pilot's 
control inputs. This matrix is a direct consequence of the weigh- 
tings on control-rate terms in Equation (4); if G = 0, t = 0. 

— — N ~ 

For single-control problems, simple, wide-band control dynamics 
and highly- responsive force-manipulators we have found that 
selecting g to yield a t n of approximately .1 sec. results in 
excellent agreement with measured data [2]. 

In more complex situations involving higher-order dynamics 
(with lower bandwidths) and less responsive manipulators, one 
might expect higher values of t N ; in a sense, the system band- 
width limitations are dominating those of the pilot. In such cases, 
it seems reasonable to select control-rate weightings on the 
same basis as the other weightings, i.e., by considering allowable 
limits for control-rates. This was the approach we took here. 

As a matter of some interest, we checked the resulting control 

* 

time constants in each case and the values obtained are pre- 
sented with the results; they lie in the range of .14 - .25 seconds, 
which is well within the range of values reported in the literature 
for neuro-motor time constants [9, 10]. 


*In a multi-control case, where the variables being controlled are 
coupled (as in the longitudinal control problem in Section 3) , the 
interpretation of control time constants is more complicated. The 
tn for a given control axis could be the corresponding diagonal 
element of the T^-matrix, or it could be an eigenvalue of T^. 

In the cases examined thus far, the eigenvalues did not differ 
significantly from the diagonal elements of Tj$. In the sequel, 
the value given for is the diagonal element of Tjj corresponding 
to the specified control. 
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(c) Pilot Randomness (Remnant ) . We assume that the 
various sources of inherent human randomness are manifested as 
errors in observing displayed outputs and in executing intended 
control movements. Thus, observation noise, v , and motor noise, 
v^, are our lumped representation of "remnant". These noises 
represent the combined effects of random perturbations in human 
response characteristics, time variations in response parameters, 
and random errors in observing displayed outputs and in generating 
control inputs . 

When the displays and controls have been optimally desig- 
ned as is the case in many laboratory situations, we would expect 
central-processing sources of pilot randomness to dominate. In 
other situations, involving more realistic displays and controls, 
display- or motor- related sources might be of greatest importance. 
All sources of remnant, however, affect the mathematical descrip- 
tions of the pilot-vehicle system in substantially the same way; 
thus, our measurements do not allow us to distinguish among 
these various noise processes. We find it convenient to reflect 
remnant largely to an equivalent observation process. Even in 
cases where the motor-noise process is needed to provide a good 
match to pilot behavior (at low frequencies) , most of the 
measured remnant is accounted for by the equivalent observation 
noise . 


For manual control situations in which the displayed sig- 
nal is large enough to negate the effects of visual resolution 
("threshold") limitations, the autocovariance of each observa- 
tion noise component appears to vary proportionally with mean- 
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squared signal level. In this situation, the autocovariance may 
be represented as 

V (t) = it P, • E{y?(t) } 
y i 

= 7T p ± Oy (t) (7) 

where P is the "noise/signal ratio" and has units of normalized 
power per rad/sec. Numerical values for of 0.01 (i.e., -20 dB) 
have been found to be typical of single-variable control situa- 
tions [ 3 , 5] . 

As noted earlier, a motor noise term is added to the pilot's 
commanded control signal to represent random errors in executing 
the intended control movements and, in addition, to account for 
the fact that the pilot may not have perfect knowledge of his 
own control activity. The motor noise is assumed to be a white 
noise, with autocovariance that scales with the control variance, 
i.e., 

= * P m. (8 > 

We have found, typically, that a value for P m of .003 (i.e., a 
"motor noise ratio" of -25 dB) yields good agreement with experi- 
mental results [2] . Throughout this study the motor-noise ratio 
was set to approximately -25 dB . 
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(d) Attention Sharing . Because the numerical value 
associated with the pilot's noise/signal ratio [P] has been found 
to be relatively invariant with respect to system dynamics and 
display characteristics, we associate this parameter with limi- 
tations in the pilot's information-processing capability. This 
forms the basis for a model for pilot attention in which the 
amount of attention paid to a particular display is reflected 
in the noise/signal ratio associated with information obtained 
from that display [4]. Specifically, the effects of attention- 
sharing are represented as 


where P^ is the noise/signal ratio associated with the i th dis- 
play when attention is shared among two or more displays, f. is 

x. l_ ^ 

the fraction of attention allocated to the i display, and P q 
is the noise/signal ratio associated with full attention to the 
display. For example, if "full attention" is represented by a 
noise/signal ratio of -20 dB, and if attention is equally divided 
between altitude-related and pitch-related information, then the 
noise/signal ratios associated with altitude- and pitch-related 
variables would be -17 dB. 

(e) Display-Related Limitations . The perceptual model 
for ideal display conditions is given by Equation (5) , with the 
covariance of the observation noise defined by Equation (7) . In 
more general situations where display characteristics are not 
ideal it is necessary to alter this model. In essence, our 
approach is to retain Equation (5) , and modify the expression 
for the observation noise covariance associated with a particular 
display variable. 
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In this study, two display limitations were important under 
certain circumstances, namely threshold limitations and the lack of 
a null-point indication. We account for these phenomena by letting 
the autocovariance for each observation noise process be 


v i (t) 


P . 
1 



\ 



,a.) 


+ 



( 10 ) 


where the subscript i refers to the l display-variable. The 
quantity K(ck, a^) in Equation (10) is the describing function 
gain associated with a threshold device 


K (a , 




where "a" is the threshold and c is the standard deviation of 
the "input" to the threshold device.* This factor is used pri- 
marily to account for threshold-type phenomena associated with 
viewing the display, but "indifference" thresholds will have an 
indistinguishable effect. Essentially, its effect is to cause 
the observation noise covariance to become greater as the 
signal becomes smaller relative to the threshold. It is worth- 
while to note that we could have introduced the threshold as a 
"gain" on the display variable y i# i.e., in the control loop 
directly. (See [6]). This is so because the optimal estimator 
then compensates for the gain and, in the process, scales the 
observation noise. Thus, it is primarily a matter of convenience 
as to whether one chooses to model threshold-related effects by 
additional elements in the control loop or by simple readjustment 
of observation noise parameters . 

♦For non-zero mean signals this expression must be modified 

(see (6] ) . 
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The term V-[ in (10) is a residual-noise covariance and 
o 

may be written as 



The residual noise term will, in many cases, be similar in 
effect to a threshold. However, it can be viewed as a separate 
parameter and used to account for observed degradation in tracking 
performance that results from lack of reference indicators [4, 7]. 

Consider the situation in which the pilot is to maintain 
the display indicator at some desired distance from a given 
reference point. The basic model for perceptual variability 
assumes that such variability is proportional to the quantity 

being estimated. Representing the mean-square as the variance 

* 

plus the square of the mean, we obtain 
V = (a 2 + y 2 ) P 

y y y 

where the mean signal value y is simply the desired distance of 
the signal from the displayed reference point. If we treat the 
quantity y as an equivalent rms "residual noise" (i.e., let 
c = y) , this expression reduces to the model for observation 
noise given earlier (with zero threshold) . 

*In most tracking studies y = 0 and we have the observation 
noise proportional to the variance of the signal. 
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Theoretical prediction of the effective residual noise 
level is somewhat more complex for the situation in which an 
indicator is to be positioned between two reference marks. A 
reasonably good approximation, however, is to let a Q be equal to 
the distance from the target point to the nearest reference. 

Model Outputs 

Once the system dynamics, displays, task requirements, 
and pilot limitations are specified it is possible to obtain 
a variety of measures of performance of both the pilot and the 
closed-loop system. These measures are obtained by solving the 
relevant optimization problems (to obtain the pilot model [2, 4, 5]) 
and by applying appropriate system analysis procedures to the 
results . 

The following measures of performance are obtained directly 
from the model or may be derived from model results: 

(a) Mean System Response . The mean time history of any 
system variable can be predicted. This is the average waveform 
that would be computed if the results of a large number of 
successive runs were averaged together at each point in space 
(time) . This ensemble average is predicted during a single run t 
however; it is not necessary to exercise the model repetitively. 

(b) RMS System Response . The rms variability about the 
mean response can be predicted for any system variable at each 
point in space. Variability arises from random inputs to the 
system such as wind turbulence and from the variability 
inherent in the pilot's response behavior (i.e., "pilot remnant"). 
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(c) Frequency-Domain Measures . Representative steady-state 
control situations can be analyzed to yield predictions of pilot 
describing functions, "open-loop" describing functions, and power 
spectral density curves for the various system variables. In 
addition, spectra for pilot remnant can be obtained. 

(d) Estimation Performance . The mean and variance of the 
estimate of the state of the system and of the error in that 
estimate (see Figure 1) . These quantities are useful in analy- 
zing the decision-making role of the pilot [12] . They are also 
useful for assessing monitoring performance with given display 
configurations . 

(e) Allocation of Attention and Workload. The model 
for attention sharing may be used to predict the effects of 

* 

task interference and to analyze workload requirements [4] . 

To predict what happens on a specific task when only partial 
attention may be paid to it, Equation (9) is used to establish 
the appropriate observation noise-signal ratio and the model 
equations are solved for this value. If the pilot's allocation 
of attention is unknown beforehand, model solutions may be used 
to determine the optimum allocation of attention, which, in line 
with the fundamental optimality hypothesis, may be taken as a 
prediction of the pilot's allocation. 

The model of task interference lends itself straight- 
forwardly to the prediction of the amount of "workload" associa- 
ted with a given task. We define the "workload index" (WI) 

*This approach may also be used as an approximate indication 
of the effects of visual scanning (see, e.g., [13]). 
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as the fraction of the controller's capacity that is required to 
perform a given task to some specified, or criterion, level of 
performance . 

Workload Index = WI = P /P (11) 

o c 

where P is the maximum noise/signal ratio that can be tolerated 
while performance is maintained within the criterion level. This 
metric can be predicted quantitatively with the existing implemen- 
tation of the optimal-control model. The procedure is identical 
to that for predicting task interference: once the model is 

"calibrated" for single-axis behavior (either by doing a simple 
experiment or by using nominal values of parameters that have been 
found to match previous data) , a curve of performance score versus 
observation noise ratio is obtained. By relating the observation 
noise ratio to fraction of capacity, a quantitative value of work- 
load may be determined. 

(f) "Window" Performance. Because the system is assumed 
to be linear and the disturbances Gaussian, all system variables 
are Gaussian. Consequently, the computation of the mean and 
covariance of all system variables provide sufficient information 
for determining the joint probability density functions of these 
variables. It is therefore possible to predict the probability 
that a variable will exceed a given value at a given time. In 
particular, the probability of a "missed approach" can be computed 
by establishing appropriate "approach windows" [14] . A single 
"run" of the optimal control model provides the probability of a 
missed approach for a given disturbance. It is also possible to 
obtain the probability of a missed approach "averaged" over all 
possible wind conditions. The procedure is described in Chapter 3. 


20 



Report No. 2484 


Bolt Beranek and Newman Inc. 


STEADY-STATE ANALYSIS OF STOL APPROACH 

In this chapter we analyze approach performance of a STOL 
aircraft using a display similar to the STOLAND EADI [7] . Steady- 
state analysis of performance at the decision height (approxima- 
tely 230m from touchdown) is the basic tool employed. Both 
longitudinal and lateral control are examined. Attentional 
allocation, window performance and pilot workload are all 
analyzed. Potential improvements of the display configuration 
are also explored. 


Vehicle Dynamics 

The vehicle considered is the Augmentor Wing Jet STOL 
Research Aircraft (AWJSRA) . The aircraft is assumed to 
be on the (7.5 deg) glide-slope at the decision height (30m), 
with a nominal airspeed of approximately 31m/s (60 knots) . 

Linearized perturbation equations for longitudinal and 
lateral-directional dynamics were used throughout. For longitu- 
dinal control unaugmented aircraft dynamics were considered. 

The pilot was assumed to control the elevator and "nozzle" 
(thrust vector) in a continuous manner, whereas the throttle 
was assumed to remain fixed at its trim setting. In the 
lateral case we considered augmented (SAS-on) dynamics. 

For this condition, we assumed that all turns were 
coordinated (i.e., zero sideslip). Longitudinal and lateral 
dynamics are presented in appropriate state-variable form in 
Appendix C. 
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Display-Related Parameters 

An abstraction of the relevant EADI display of status 
information is shown in Figure 2. This display provides the 
pilot with glide path and localizer errors as well as attitude 
information. From such a display the pilot can also obtain the 
rates of change of these variables. Although we have not shown 
an airspeed error indicator in Figure 2, the pilot is displayed 
this quantity and we will assume that airspeed error is available 
in our analysis. 

Effective visual thresholds were computed for the air- 
craft at the 30-meter decision height. On the basis of previous 
analysis of approach performance [6] , an "indifference threshold" 
of 0.1 degrees visual arc was associated with perception of 
height error. Previous analysis of pilot remnant data [5] 
suggested thresholds of 0.05 degrees visual arc for other 
indicator displacements and 0.18 arc-degrees/second for indi- 
cator-rate quantities. Display gains given in [7] were used to 
convert thresholds into units related to system quantities. 

Non-zero rms residual noise terms were associated with 
height and sink-rate information. Since the status display shows 
the height window of +3.7 meters, a residual noise on height per- 
ception was needed to account for the lack of an explicit zero 
reference. To approximate the effects of the non-zero reference, 
we simply set the value of the residual noise equal to the amount 
of the reference offset (i.e., 3.7 meters). The residual noise 
on sink-rate information was included to account, in a rough way, 
for the resolution limitations of the SBILS (Scanning Beam 
Instrument Landing System) . As we show later, the latter 
noise term had no appreciable effect on predicted performance. 


22 



Report No. 2484 


Bolt Beranek and Newman Inc. 


( ^ 

GLIDE SLOPE AND 



FIGURE 2. Display of Status Information 


23 



Report No. 2484 


Bolt Beranek and Newman Inc. 


A non-zero rms residual noise term was also associated with 
lateral offset error. This was set to a value of 5.4 meters, 
corresponding to the lateral dimensions of the "window". 

Thresholds and residual noises for all displayed vari- 
ables are summarized in Table 1. 

Di s turbances 


Zero-mean, random gusts provided the disturbances for 

this analysis. These were generated by passing white-noise 

through first-order filters, such that the gust spectra 

approximated those of Reference 16, for the 30m altitude. 

Specifically, the three components of turbulence (u , v , w ) 

* 9 y y 

were generated by 


g 


+ 
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g 
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(12a) 

(12b) 


(12c) 


where 

= L v = 205m ; = 30.5m 



( 12d) 
(12e) 


♦Effective angular rates due to gusts 





were neglected. 
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Table 1 

DISPLAY- RELATED PARAMETERS 


Variable 

Threshold 

RMS Residual 
Noise 

h (m) 

0.48 

3.7 

• 

h (m/s) 

0.85 

1.1 

6 (deg) 

0.22 

0.0 

q (deg/s) 

0.78 

0.0 

u . (m/s) 

0.14 

0.0 

y (m) 

.28 

5.4 

• 

y (m/s) 

. 52 

0 

★ 

ci) (deg) 

0 

0 

. * 

| <+, (deg/s) 

o 

0 

' 





♦Studies of longitudinal control showed that the thresholds on 
attitude could be ignored, so in the subseguent analysis of 
lateral performance the thresholds associated with viewing the 
bank angle indicator were set to zero. 
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and V : 31 m/s is the aircraft velocity. The longitudinal gust 
spectrum corresponding to (12a) is identical to that in [16] , 
whereas the vertical and lateral gusts spectra are approxima- 
tions to those of [16] . 

Gust intensities were assumed to be distributed according 
to a Rayleigh distribution [16]. For most of the analyses, gust 
intensities were either 

a =2.1 m/s, a = .82 m/s, a v =2.1 m/s 

u g w g g 

or 

a = .8 m/s, a = .31 m/s, o = .8 m/s 

U W V 

g g g 

According to [16], the probability of encountering gust amplitudes 
of the magnitude of (13a) , given that any turbulence is encountered 
is about 1%, whereas the magnitudes of the level of (13b) have 
about a 50% probability of being encountered. Thus, (13a) corres- 
ponds to a "worst-case" disturbance and (13b) to a "median" dis- 
turbance . 


(13a) 

(13b) 


Task Requirements 

As noted in Chapter 2, task reguirements are stated in terms 
of cost weightings corresponding to "allowable" limits of variables 
The "limits" and weightinas corresponding to the longitudinal 
and lateral control tasks are given in Tables 2A and 2B, res- 
pectively. Height, airspeed and lateral deviation limits were 
based on Category II "window" specifications, control and control- 
rates were determined largely from physical considerations, and 
the remaining weightings were based on human preferences. Also 
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Table 2 

LIMITS AND COST FUNCTIONAL WEIGHTINGS 
Longitudinal 


Variable 


"Limit" 

3.7 

(m) 

1.1 

(m/s) 

6.0 

(deg) 

2.6 

(m/s) 

9.0 

(deg) 

29. 

(deg) 

50. 

(deg/s) 

100. 

i 

(deg/s) 


Weiahtina 


0.07 3 
0.83 
0.028 
0.0 
0.15 
0.012 
0.0012 
0.0004 
0.0001 


Lateral 


Variable 


"Limit" 


5.4 (m) 

20. (deg) 

30. (deg) 
60. (deg/s) 


Weighting 


.034 

0 


0 

.0025 

.0003 
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shown in Tables 2A and 2B are the t^'s (i.e., the corresponding 
diagonal elements in the T^ matrix) that result from the selected 
weightings . 

Pilot N o i s e / S i g n a 1 Ratios 

In all the analyses performed here the motor noise/signal 
ratio for each control variable was approximately -25 dB. This 
value was based on previous results [2] . The "base" observation 
noise/signal ratio (P ) was treated as a measure of overall 
attention to the task being considered and was made the same 
for all display Quantities in the task. The noise/signal 
ratio associated with a given display variable was obtained by 

dividing this base level by the corresponding fraction of atten- 

* 

tion. The value for P q will be a parameter in our studies of 
workload. On the basis of previous analysis of remnant data 
[2, 3], an observation noise/signal ratio of -20 dB was consi- 
dered as the "nominal" value. We will discuss this choice 
further in connection with workload analysis. 

Analysis of Longitudinal Performance 

The optimal-control model was used to analyze performance 
in longitudinal control with the dynamics of Appendix C, the 
display parameters of Table 1 and a cost functional corresponding 
to the "weightings" in Table 2A. 

*The "total" nois"e associated with a display variable is then ob- 
tained by modifying this level to account for threshold and 
residual-noise effects as in Equation 10. 
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Allocation of Attention . — A "sensitivity" study was 
performed to determine the "optimal" allocation of attention 
among the three longitudinal display variables, heiaht error 
(h) , pitch (6), and airspeed (u) . 

Table 3 shows the fractional allocation of attention to 

* 

these display variables. This distribution of attention was 
found to yield the lowest total predicted cost, subject to 
the constraint that the fractional attentions sum to unity. 

It is important to note that performance was relatively in- 
sensitive to attentional allocation, and one could vary the 
allocation over a fairly wide range without affecting the pre- 
dicted results appreciably. The values for attention given in 
Table 3 were used for the remainder of the longitudinal analysis 

Table 3 

ALLOCATION OF ATTENTION AMONG LONGITUDINAL DISPLAYS 


Variable | 

, 

_ - L 

h , h 

CD 

-Q 

l 

u ! 

r 

i 

Attention 

.35 

.1 

1 

. 55 


I 


*The pilot is assumed to obtain derivative information from a 
given display indicator without additional attention being 
recruired. Thus, the fractional attentions to height-related 
Quantities as a group, pi tch- related quantities as a group, 
and airspeed sum to unity. 
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Per f orma nee . — Table 4 shows the predicted rms perfor- 

mance scores for the various display and control variables for 
the "worst-case" wind condition, with attention allocated "optimally" 
and all other model parameters set at "nominal" values. Also 
shown are the "limits" associated with each variable (repeated 
here for convenience) and the probability that the magnitude of 
a given variable will exceed its limit. The latter was computed 
as twice the integral of the Gaussian probability density function 
from the "limit" to infinity. (The factor of 2 was included to 
account for both positive and negative overvalues.) 

Four variables have a probability greater than 0.001 of 
exceeding their respective limits: height (11%), sinkrate (10%), 

airspeed (5%), and nozzle (2%). If we assume that the probabi- 
lity distributions for height and airspeed errors are essentially 
independent, we can use the marginal probabilities to compute a 
joint probability of 16% for simultaneously exceeding the height 
and airspeed windows. 

Table 4 

PERFORMANCE SCORES FOR THE "NOMINAL" CONTROL SITUATION 

(0 = 2.1 m/s, a = .82 m/s, t = .2 sec., P = -20 dB) 

u w o 


Variable 

"Limit" 

RMS Error 

P ( | Error | > Limit) 

h (m) 

3.7 

2.3 

0.11 

h (m/s) 

1.1 

.68 

0.10 

0 (deg) 

6.0 

1.8 

0.001 

q (deg/s) 


1.2 

0.0 

Uj(m/s) 

2.6 

1.3 

0.05 

6^deg) 

9.0 

2.1 

<0.001 

j <S N (aeg) 

29 . 

7.4 

<0.001 

| 6 e (deg/s) 

50. 

12. 

0.02 

) 8 (deg/s) 

N 

100. 

21. 

<0.001 


30 


V 
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The "nominal" control situation described above corresponds 
to the worst-case (1%) wind condition and one can generally expect 
better performance than indicated in Table 4. A more reliable 
indication of system performance would consider the "range" of 
possible wind conditions. In order to obtain such a measure of 
pilot/vehicle performance, it is necessary to specify how performance 
degrades as a function of oust intensity. If pilot-related noise 
processes were not considered, rms errors would simply vary 
linearly with rms gust amplitudes in the absence of other external 
inputs to the system. However, because thresholds and resolution 
limitations are included in the model, a more complex relationship 
can be generally expected. 

The model was used to predict rms performance measures for 
five sets of gust intensities, ranging from the 80%-wind to the 

0.25%-wind. The effects of gust intensity on rms height, sink 

* 

rate, and airspeed errors are shown in Figure 3. All three of 
the perf ormance-versus- intensity curves are very nearly linear 
and straight-line approximations are obtained by connecting the 
rms levels corresponding to the 50%~and 1%-winds. Because of 
the consideration of threshold effects, the linear approximations 
for height and sink rate errors have positive zero intercepts. 

We have shown thus far how to predict svstem performance 
as a function of rms gust intensity. Since all svstem variables 
are assumed to be Gaussian, we can use the predicted rms perfor- 
mance scores to obtain the probability density function for any 
variable, conditioned on the rms gust level. It is worth noting 
that although rms performance scores vary linearly with gust 
intensity, the probabilities of exceeding approach window limits 
do not (Figure 4). 

^"Note t h at" sp i c if y i ng — a also specifies a via Eguation (12e) . 

u g w g 
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RMS LONGITUDINAL WIND GUST (METERS/SEC) 


FIGURE 3. Effect of Gust Intensity on Performance Scores 




IT* 



0.4 0.8 1.2 1.6 2.0 2.4 

RMS LONGITUDINAL GUST INTENSITY, <r , METERS/SEC 

u g 


FIGURE 4. Effect of Gust Intensity on "Window" Probabilities 
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In order to predict a measure of system performance that 
is averaged over all possible wind conditions, the joint pro- 
bability distributions of system variables and rms wind inten- 
sities has to be integrated. Thus, if one wishes to predict 
the overall probability of a height error greater than the 
maximum allowable value, the following expression must be 
evaluated : 


P ( | h 


h L ) = 2 






oc 

- OO 

= 2 • P / p (a ) 

g J u g 

(p(h|a u )dh 

0 

A t 9 



(14) 


where p(a Ug ) is the probability density function of the rms 

longitudinal gust, p(h|a u ) is the conditional Gaussian pro- 

g 

bability density function for height error, and is the 
probability of encountering clear-air turbulence.* The factor 
of 2 is included so that negative as well as positive over- 
values will be accounted for. 


The density function p(a u _) may be taken as that of a 
Rayleigh distribution [16]. 


P(°u ) 
g 



-1/2 c 
e 



(15) 


♦According to [ 16] , there is a probability of about 0.2 that no 
turbulence will be encountered at an altitude of 30 meters. 
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where c = .703 m/s. If rms-height (or any other variable of 

interest) was directly proportional to a u , i.e., if 

y 

a^/a u = K = constant (16) 

g 

then Equation (14) could be evaluated analytically.* But, because 
of thresholds, Equation (16) does not hold. Conseouently , a numeri- 
cal integration technique was employed to evaluate Equation (14) 
and a similar expression for airspeed. Linear approximations to 
the curves of rms performance vs. rms gust intensity were used. 

For the case considered above (i.e., P Q =-20dB) , the overall 
probability of beinq out the height-airspeed window turned out to 
be about 1.2%. 

Workload Requirements . — Pilot workload requirements 
under the various gust conditions can be determined from the 
relationship between performance and noise/signal ratio. Re- 
call the model for pilot attention, in which 

f = P /P 
o' 

where f is the fraction of attentional capacity devoted to the 
task, P is the noise/signal ratio achieved by the pilot when 

*The easiest way is to inteqrate by parts. Let s = a /c , 

u g 

(14) may be written as 

CO 

/ [u]dv, which can be 

o 


t = h/o^ an< ^ H = h L /Rc. 


-t2/2 


°o r °o, ^ 

"o 1 H/s /2 tT ^ _ H 

evaluated to yield ie 


dt 


Then Eauation 

1 - S 2/2, 


se 


’ds = 
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performing the task, and P Q is a reference ratio. Building on 
this model for attention, we have defined a "workload index" as 
the amount of attention required to achieve a specified criterion 
level of performance on the control task. Thus, 

Workload Index = P /P 

o c 

where P c is the maximum noise/signal ratio that can be tolerated 
while performance is maintained within the criterion level. 

In order to predict the workload index, we have to specify 
some relevant performance measure, the required level of perfor- 
mance, and the "reference" noise/signal ratio P . Ideally, we 
would like P Q to correspond to full attention. Unfortunately, 
we cannot conduct an experiment in which the pilot is guaranteed 
to use his total information-processing capability. Therefore, 
we let P Q correspond to the noise/signal ratio (namely, .01 or 
-20dB) obtained in a standardized laboratory situation in which 
the pilot is motivated to minimize his tracking errors. We know 
that this value does not correspond to "full capacity", because 
significantly lower noise-ratios have been found experimentally 
[4], However, based on our laboratory experience, P Q = -20dB 
does appear to correspond to a high workload condition, and 
"operation" at this level for any prolonged time would un- 
doubtedly be unacceptable. Of course, when we are interested 
primarily in the relative change in workload requirements from 
one situation to the next, the value for P Q is not too critical; 
this is the case in this study. 
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We consider the probability of a missed approach as the 
relevant measure of system performance for assessing workload. 

We define a "missed approach" as a situation in which either 
the height or the airspeed error exceeds the maximum allowable 
value as defined by the "limits" shown in Table 2A. The 
specification of a desired level of performance is somewhat 
arbitrary, but we suspect that a reasonable criterion lies 
between 1% and 5% (i.e., a probability of a successful approach 
between 95% and 99%) . 

The relation between noise/signal ratio and the pro- 
bability of exceeding the height-airspeed window is shown in 
Figure 5 for the 1% and 50% wind conditions. Also shown are 
the curves for exceeding the height and airspeed limits indi- 
vidually for the 1% wind condition. The allocation of attention 
among the height, pitch and airspeed display variables was held 
at the levels shown in Table 3 as the overall noise/signal ratio 
was changed. Because of the low sensitivity of performance to 
attentional allocation, we do not expect that significantly 
different performance scores would have been predicted had we 
re-optimized the allocation at each stage in the analysis. 

For the 1% wind condition, the predicted probability of 
a missed approach is about 10% for a noise/signal ratio of -26dB 
( f a;4 . 0 ) — the lowest ratio explored. Thus, the workload index 
will be well above 4.0 for a criterion level in the range of 1% 
to 5%. Since noise/signal ratios as low as -26dB have been 
found experimentally only in extremely demanding control situa- 
tions (unstable dynamics) [4] , we would conclude that the work- 
load requirements for the approach task in a 1% wind environment 
are excessive. 
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If the wind intensity is reduced to the 50% level, however, 
the task workload is reduced considerably. For a 1% criterion, 
the workload index is slightly less than 0.5, which represents 
an achievable (if not comfortable) level of attention. If we 
relax the criterion to 5%, the workload index is reduced to 
about 0.23. Thus, for a 5% criterion probability, the work- 
load requirements for the task range from low-to-moderate for 
the 50% wind condition to excessive (if not impossible) for 
the 1% wind condition. 

The contributions of height and airspeed to the overall 
probability of a missed approach can be seen from the corres- 
ponding curves shown in Figure 5. If the pilot were to operate 
at a noise/signal ratio of -26dB, height and airspeed errors 
would contribute about equally to the 10% probability of a 
missed approach in the 1% gust environment. Since height regu- 
lation is considerably more sensitive to pilot-induced noise 
than is airspeed regulation, however, the relative contribution 
of height errors increases markedly as attention is relaxed. 

For the 50% wind condition, a missed approach is caused by 
height errors only. For the range of noise/signal ratio 
explored, the probability of excessive airspeed errors is 
less than 0.1 percent. 

Figure 6 shows the predicted probability of a missed 
approach as a function of observation noise/signal ratio, avera- 
ged over all possible gust intensities. If the pilot operates 
consistently at a noise/signal ratio of -20dB, he will achieve 
an overall probability of a missed approach of about 1.2 percent. 
Thus, the average workload index is about 1 if we specify a 
criterion performance level of 1%. If a probable missed approach 
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of 5% is adopted as the performance criterion, the average work- 
load index is reduced to about 0.25. Figure 7 shows the trade- 
off between workload and performance (probability of a missed 
approach) for the longitudinal task and the given display. 

Effects of Display Parameters . — The foregoing analysis 
was conducted for the nominal display condition (see Table 1) in 
which visual thresholds, other types of display-related noise, and 
attention-sharing requirements served to limit performance. In 
order to explore the benefits that could be obtained by reducing 
various display-related limitations, model analysis was performed 
in which these limitations were removed in stages. The display 
configuration was adjusted as follows: 

Condition A: nominal display condition. 

Condition B: removal of the residual noise term 

associated with height perception. This condition 
corresponds to providing the pilot with an explicit 
zero reference for glide-slope error. 

Condition C : removal of the residual noise term on 

sinkrate to correspond to a noise-free SBILS pre- 
sentation. 

Condition D: zero threshold for sinkrate perception. 

Condition E: zero thresholds on all remaining display 

variables . 

Condition F: No modification of noise/signal ratios to 

account for attention-sharing. This condition corres- 
ponds to an idealized integrated display. 
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The model conditions described above were cumulative. 
Once a deviation from the nominal condition was made, it re- 
mained in effect for the remaining parameter adjustments. 

Not all of the above conditions are intended to be 
physically realizable. For example, a zero visual threshold 
corresponds to infinite display gain. Thus, this analysis 
serves to indicate the upper bounds on performance improvement 
that can be expected from modifying various display-related 
variables . 

Except for changes in display-related parameters, the 
model parameters were selected for the nominal condition. Only 
the 1% wind condition was considered. The results of this 
analysis are presented in Figure 8. For each display condition 
defined above, probabilities of excess excursions were pre- 
dicted for height, sinkrate, and airspeed. In addition, the 
probability of a missed approach (the combined probability of 
exceeding either the height or airspeed limits) was computed. 

Performance along all dimensions except airspeed regula- 
tion improves considerably as the display- related limitations 
are removed. The greatest improvement is predicted when the 
effective visual threshold on sinkrate information is set to 
zero. Small-to-moderate benefits accrue from simulating a 
zero reference for glide slope error and from neglecting the 
effects of attention-sharing. Negligible improvements result 
when remaining threshold terms are set to zero. 

Overall, the predicted probability of a missed approach 
decreases from about 15% for the nominal condition to about 
6% for the idealized display condition. Thus, predicted work- 
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load remains relatively high even for the idealized display condi- 
tion. (Results were obtained for a noise/signal ratio of -20dB.) 
It must be remembered, however, that these predictions correspond 
to the 1% wind condition only, and that meaningful improvement 
might be observed if we consider performance averaged over 
all wind conditions. 

It is not surprising that performance should improve con- 
siderably when the threshold on sinkrate is removed. Comparing 
Tables 1 and 4 we observe that the predicted rms sinkrate error 
is less than the assumed threshold. Thus, the pilot can be 
expected to obtain no useful sinkrate information from the 
velocity of the glide-slope error indicator for the nominal 
display configuration. Since sinkrate information is useful in 
predicting future flight-path errors, performance should be 
improved if the effective threshold on sinkrate can be suffi- 
ciently decreased. (Note that if we attempt to accomplish this 
by increasing the gain of the glide-slope indicator, we run the 
risk of unacceptably large excursions in the high-gust situation.) 
Experimental results recently obtained at Ames also support the 
prediction that a decreased threshold on sinkrate will lead to 
better height regulation. In a study of approach performance 
with a pictorial display, performance was found to improve con- 
siderably when the pilot was shown a predicted touchdown point 
[6] . In effect, a presentation of the aim point error provided 

the pilot with a high-gain display of sinkrate error. 

/ 

The source of the improvement in performance with the 
idealized display is, more generally, an improvement in moni- 
toring performance. The standard deviation of the estimation 
error predicted by the model normalized with respect to the 
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predicted standard deviation of the signal gives an indication of 
how well a pilot may be expected to monitor a given variable. 

The results for the status and the idealized display are presen- 
ted in Table 5. It can be seen that monitoring of all variables 
is improved. The improvement is greatest for height errors and 
least for speed errors. The relatively small improvement in 
monitoring airspeed is reflected in the lack of substantial 
improvement in airspeed regulation. 

Table 5 


MONITORING PERFORMANCE FOR STATUS 
AND IDEALIZED DISPLAYS 

(1% wind, P Q = -20 dB) 


Variable 

0 

S 

tatus 

CT e 

a e/o 

a 

Idealized 

°e 

° e/a 

h (m) 

2.31 


.96 

.415 

1.53 

.17 

.111 

h (m/s) 

.68 


.41 

.603 

.55 

.18 

. 327 

6 (deg) 

1.81 


.47 

.261 

1.77 

.22 

.124 

q (deg/s) 

1.22 


.20 

.164 

1.91 

.17 

.089 

u(m/s) 

1.34 


.78 

.586 

1.30 

.70 

. — — ■ - .. < 

.537 


a = Standard Deviation of Signal 

a = Standard Deviation of Estimation Error 
e 
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In summary, the results of the analysis of display-related 
parameters shows that significant improvement in system perfor- 
mance can be expected if the pilot is presented with a zero 
reference for glide-slope error, if the ability to obtain sink- 
rate information is considerably enhanced, and if the effects 
of attention-sharing can be avoided. One way to approach this 
ideal situation is by providing the pilot with a flight director. 
With this display the pilot need look at only one or two display 
indicators, zero references are provided, and the gains on the 
various display quantities are optimized so that the pilot is 
shown the appropriate "mix" of inputs on a display that is 
suitably scaled. Thus, the model predictions obtained under 
Condition F (described above) give a good indication of the maxi- 
mum benefits than can be obtained from a flight director. Accor- 
dingly, the results of this analysis will be used later in Chapter 
4 to evaluate proposed flight director laws. 

Analysis of Lateral Performance 

Lateral performance using the status display described 
earlier was also analyzed. Vehicle dynamics for the analysis are 
tabulated in Appendix C and correspond to a simplified situation 
in which all turns are assumed to be coordinated; lateral gusts 
are as described above; display parameters and cost weightings 
are given in Tables 1 and 2B above. 

Stability Augmentation . — A brief preliminary investiga- 
tion of SAS-off lateral dynamics (also tabulated in Appendix C) 
was conducted for the basic display configuration. Figure 9 
compares lateral performance vs. observation noise/signal ratio 
(attention) for the SAS-on and SAS-off cases. It may be seen 
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FIGURE 9. Effect of Stability Augmentation and Attention on 
Lateral Control Performance 
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that at low noise/signal ratios it is actually possible to achieve 
smaller lateral errors with the SAS-off. However, the SAS-off 

condition is extremely sensitive to observation noise/signal 

★ 

ratio and is, consequently, auite demanding. It is auite 
likely that this sensitivity is what makes the SAS-off confi- 
guration unacceptable. In the remainder of our analysis we 
consider only SAS-on dynamics. 

Allocation of Attention . — A sensitivity study was per- 
formed to determine the optimal allocation of attention between 
localizer and bank angle displays. The results, for two values 
of P , are presented in Table 6. At the lower noise/signal 
ratio the "optimal" allocation is about 70-75% attention to 
the localizer whereas at the higher noise/signal ratio it appears 
that full attention should be devoted to the localizer. However, 
in both cases the results are auite insensitive to allocation of 
attention. In the remainder of the analysis we assumed that 75% 
of the pilot's attention would be devoted to the localizer 
display. 


Performance . — Lateral tracking performance was computed 

for the worst-case (1%) and median (50%) winds. Scores for the 
1%-wind and the nominal observation noise/sianal ratio of -20dB 
are shown in Table 7, along with the "limits" used in selecting 
the cost weightings and the probability of exceeding these 
limits . 


The (assumed) linear relation between scores and gust 
intensity was determined by passing a straight line through the 
results for the worst-case and median winds. This line was then 
used in Equation (14) to compute the probability of exceeding the 

"*In conducting the computer analysis it was also found to 
be highly sensitive to motor noise. 
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Table 6 

EFFECTS ON LATERAL PERFORMANCE OF SHARING ATTENTION 
BETWEEN LOCALIZER AND BANK ANGLE INDICATOR 


Attention 

to 

Localizer 

P = 
o 

-20dB 

0 

ii 

-14dB 

Cost(J) 

c y (m) 

Cost (J) 

Oy (m) 

.25 

1.57 

5.70 

- 

- 

.5 

1.37 

5.22 

2.43 

7.17 

.6 

1.35 

5.15 

- 

- 

.65 

1.34 

5.11 

- 

- 

.7 

1.33 

5.08 

2.21 

6.75 

in 

r- 

• 

1.34 

5.09 

2.18 

6.68 

.8 

1.34 

5.09 

2.16 

6.62 

.85 



2.14 

6.58 

.9 


- 

2.14 

6.56 

.95 


- 

2.13 

6.54 

1.0 

- - 

1.41 



5.17 i 

f 

i 

1 
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Table 7 

LATERAL PERFORMANCE SCORES 

<c v = 2.1 m/s) 

g 


Variable 

Score 

"Limit” 

Prob. of Exceedance 

a y (m) 

5.08 

5.4 

.288 

°y (m/s) 

1.99 

- 

- 

°4> (deg ) 

5.06 

20 

.000 

(deg/s) 

3.28 

- 

- 

a x (deg) 
°w 

10.5 

30 

.004 

Of (deg/s) 
w 

19.9 

60 

.003 
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lateral window (5.4m) averaged over "all" winds. This was done 
for several levels of total attention devoted to the lateral 
task. The results for the "overall" average performance, as 
well as those for the 1% and 50% winds are presented in 
Figure 10. 

Comparison with longitudinal results (Figures 5 and 6) 
reveals that the lateral control task, even with the SAS-on, is 
more difficult. The probability of missing the lateral window 
when averaged across all winds is 1.5 - 3 times as great as that 
for missing the longitudinal window at all levels of attention 
investigated. In other terms, if we establish a 95% probability 
of a successful lateral approach as a criterion level, the 
Workload Index for the lateral task, for the 50% wind is about 
.4, and it is about .7 when the average of all winds are consi- 
dered. For the 1% wind, it does not appear that a success 
probability of 95% is achievable within the limits of behavior 
that we have observed heretofore. 


Display Related Parameters . — An analysis of the sensi- 
tivity of performance to changes in display parameters, completely 
analogous to that for longitudinal control, was conducted. The 
display configuration was adjusted as follows: 


Condition A: 


nominal configuration. 


Condition B: 


removal of residual noise associated with 
lateral error (providing a zero-reference) . 


Condition C: 


zero threshold for lateral error-rate. 
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Condition D: zero threshold for lateral error. 

Condition E: no modification of noise/signal ratios 

for attention-sharing (display integ- 
ration) . 

The "worst-case" wind condition was used and the para- 
meters relating to the pilot were "nominal". Relevant scores 
and monitoring performance (as defined earlier) are given in 
Table 8. 

Table 8 

EFFECT OF DISPLAY PARAMETERS ON LATERAL PERFORMANCE 

( v =2.1 m/s, P_ = -20dB) 
g o 


Condition 

A 

B 

C 

D 

E 

o y (m) 

5.09 

5.05 

4.89 

4.87 

4.59 

°y (m/s) 

2.0 

1.99 

1.96 

1.95 

1.88 

°<f> (deg) 

5.04 

5.03 

4.97 

4.96 

4.71 

o r (deg) 

10.6 

10.6 

10.4 

10.4 

9.8 

°w 






a /a 
Ye Y 

.19 

.16 

.14 

.14 

.12 

a • /o’ 

.45 

.45 

.42 

.42 

.41 

y e y 






o, /a, 
<f> e 4> 

.12 

.12 

.12 

.12 

.10 


Performance improvements with display changes are not 
too dramatic with the cumulative improvement in tracking perfor- 
mance being about 10%. Of the changes made, only two had any 
significant effect; namely, removing the lateral error rate 
threshold and removing the necessity for attention sharing, 
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although zeroing the residual noise on error did result in a 
significant improvement (15%) in monitoring that variable. This 
suggests that performance improvements with a lateral flight 
director might not be substantial. It should be remembered, 
however, that this analysis was conducted for a high-wind, 
high-workload (P = -20dB) condition. The advantages of a 
flight director may prove to be most significant in situations 
where less "attention" is involved. Alternatively, it might 
allow achievement of similar performance at reduced workload, 
i.e., it might reduce the workload index. 

Combined Longitudinal and Lateral Tasks 

Thus far, we have considered longitudinal and lateral 
tasks separately. This is possible because of the decoupling 
inherent in the linearized perturbation equations. Of course, 
the pilot must share his capacity between the longitudinal and 
lateral tasks, which implies some interference and a degra- 
dation in performance on each task. This interference may be 
treated within the same framework of the model of attention 
presented earlier. Moreover, one can make direct use of the 
results obtained in the separate studies of longitudinal and 
lateral control. 

To account for the interference, we define a combined 
cost functional 

J TOT = J LONG + J LAT 

where and J T . m are the cost functionals for the lonaitudinal 

LONG LAT 

and lateral cases, respectively (see Table 2) . The combined cost 
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functional is meaningful because of the manner in which the 
separate cost functionals were defined and normalized. Now, 
if fLQNG is t ^ ie fraction of attention devoted to the longitu- 
dinal task, then the fraction devoted to the lateral task is 

f LAT = 1 ^ LONG 

Using the results of the earlier "attention" studies, it is 
possible to determine how attention should be shared between 
the two control modes so as to minimize ^ T q TAL “ (We assume 
that the attention sharing among displays within a given mode 
is unaltered.) 

Figure 11a shows how performance varies with attention 
to the longitudinal task for the 1%-wind and Figure lib presents 
the same data for the 50% wind. These results assume that 
P Q = -20dB corresponds to full attention. It may be seen that 
attention should be divided nearly eaually between the two tasks 
with a slight bias to the lateral task in the high-wind case. 

Figure 12 shows the effects of attention-sharing on 
window performance, averaged over all winds. Using the probability 
of a missed approach as the measure of performance leads to the 
conclusion that approximately a 40/60 split of attention between 
longitudinal and lateral tasks is optimal, thus confirming the 
results obtained using J T0T . The corresponding overall probabi- 
lity of a missed approach (i.e., a miss on height or airspeed 
or lateral position) is about 8%. 

Finally, if we assume that the pilot shares his capacity 
equally between the two tasks (which is not far from the optimal 
strategy) , we can use the earlier results to obtain missed approach 
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FIGURE 11a. Effect on Optimal "Costs" of Attention Sharing 

Between Longitudinal and Lateral Tasks (1%-Wind) 
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FIGURE 12. Effect on Approach Performance of Attention Sharing 
Between Longitudinal and Lateral Tasks 
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probabilities as functions of the relative attention devoted to 
the tracking task as a whole (assuming, for convenience, "full" 
attention r -20dB, as before).* The result is plotted in 
Figure 13. This figure emphasizes the difficulty of the task. 
When all winds are considered, it does not appear possible to 
achieve a 95% approach success probability, at least within the 
range of pilot workload that is assumed acceptable. Even for 
the 50% wind condition, a success probability of 95% implies a 
Workload Index of about .9, hardly a desirable situation. 

V SUMMARY 


The performance of the AWJSRA with an "unaugmented"EADI-display 
of status information was analyzed with the steady-state optimal 
control model for pilot-vehicle analysis. Both "window" per- 
formance and pilot workload were investigated for a range of 
turbulence conditions. 

The results indicate that with the basic display the 
overall task is cruite difficult. When the median wind level is 
considered, a 95% success probability for approach requires a 
high workload. If performance is averaged over all possible 
winds, such a success-probability does not appear to be at- 
tainable within a reasonable range of workload. 


The lateral-directional task seems to be considerably 
more difficult than the longitudinal control task, even though 


*Thus, in this analysis, the "base" observation noise/signal ratio 
for the longitudinal and lateral "subtasks" is -17dB. This, in 
turn, implies that we can't simply add the "relative" attentions 
for the longitudinal and lateral tasks, to obtain the relative 
attention for the combined task (that is, of course, the result 
of interference) . It should also be mentioned that combined 
performance may be computed for any split of attention. The 
"equal-split" allows us to use earlier computations directly, 
i.e., without interpolation. 
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stability augmentation is provided for lateral control. For a 
95% probability of being within the respective approach window, 
the lateral task has a workload index of about .63, as compared 
to a workload index of about .25 for longitudinal control. 

Potential improvements to the basic status display were also 
explored. The greatest effects were observed when better error- 
rate (sink-rate, lateral error-rate) information was assumed, 
as might be provided, for example, by a display of longitudinal 
and lateral flight path angles. Significant effects were also 
observed when the recruirements for attention-sharing were removed. 
These improvements, as well as a reduction in pilot workload, may 
be realizable with a flight-director display. 


62 



Report No. 2484 


Bolt Beranek and Newman Inc. 


ANALYSIS OF PROPOSED FLIGHT DIRECTOR LAWS 

In the preceding chapter it was shown that the unaugmen- 
ted status display was probably unacceptable both from a closed- 
loop system performance and a pilot workload standpoint. Approp- 
riately designed flight directors might serve to alleviate many 
of the problems. In this chapter we analyze potential improve- 
ments with a set of proposed "interim flight directors" [17] 

(or, more exactly, approximation thereto) . Our analysis here is 
confined to steady-state gust regulation; it should be emphasized 
that the flight director laws we are examining were based on 
broader concerns, such as wind shear compensation. 

Longitudinal Director System 

The longitudinal director laws considered herein are 

simplifications of the interim director laws [17] , but should 

* 

retain the important features of those laws. The simplifications 
serve to reduce the computational requirements. A block diag- 
ram of the approximate implementation is given in Figure 14. 

The corresponding gains are those of [17] and are listed in 
Table 9. A discussion of the basic design rationale for this 
director system may be found in [18]. 


*Private communication with R. Klein. 
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Table 9 

LONGITUDINAL FLIGHT-DIRECTOR GAINS 

Kg = 1.0 units FD/RAD; T g = 1.0 sec 

K = 1.0 units FD/RAD/SEC ; T 1 = 3.0 sec 

K u = JO 328 units FD/METER/SEC ; T 2 = 1.0 sec 

K* = .0328 units FD/METER/SEC 

d 

s 

K D = 1 degrees/unit 

s 

K , = 1.62 units FD/METER 

d 

K* = 3.28 units FD/METER/SEC 

d 

n 

K x = 30.8 units FD/RAD-NOZZLE DEFLECTION; T 7 = 1.0 sec 

6 N 

Model Parameters 

Thresholds for the stick director (FD g ) were set at the 
basic visual threshold values obtained from remnant studies [4] ; 
given the scaling, the thresholds associated with the nozzle 
director were neglected. Inasmuch as both directors have zero- 
references, residual noises were set to zero. Pilot-centered 
parameters, of time-delay and motor-noise/signal ratio were set to 
the nominal values of .2 sec and -25dB, respectively; the "base" 
observation noise/signal ratio, P , was -20dB. 
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Task Requi rements 

The choice of a suitable cost-functional for analysis 
of the flight director displays is not so straightforward as in 
the no-director case. To choose this cost functional one must, 
in some measure, determine how the pilot interprets the task 
when a flight director is used. This will, of course, depend 
on his instructions, and his understanding and training. Four 
modes of operation suggest themselves as potential pilot 
strategies : 

i) The task is viewed as a simple compensatory tracking 
task (2-axis, here) with the pilot attempting to 
minimize mean -squared flight director error. 

ii) The task is viewed as being the same as for the 
no-director situation and the purpose of the 
director is, essentially, to provide better 
information. 

iii) The director signal is interpreted in terms of 
its component inputs and its maximum allowable 
excursion is balanced against "limits" on control 
and control-rate. 

iv) Maintain director errors small while keeping all 
other variables within bounds. 

If the first strategy is assumed, the cost weightings are chosen 
on the same basis as for laboratory tracking tasks; i.e., only 
error and control-rate contribute to the score and the control- 
rate weighting is chosen to yield a desired t n . This assumption 
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has little to recommend it but simplicity* because it ignores 
flight experience. A few trial runs were made with such a 
cost functional; they resulted in extremely small tracking 
errors but unacceptably large control inputs. 

If we assume the second strategy to be operative, then 
we use the same cost functional that was used before. This 
assumption has within it the implicit notion that the pilot will 
also use his status displays to assure that untoward vehicle 
motions are not induced by blind following of the director 
signal. Choosing this cost functional also has the advantage 
of providing for a direct, one-to-one, comparison between 
director and no-director cases. The principal objections to 
this approach would appear to be that it doesn't emphasize the 
zero-reader aspects of the use of the directors and it is less 
acceptable if one considers only director signals to be available. 

To implement the third strategy we consider the inputs 
to the flight director. For example, the stick director signal 
is scaled so that one degree of pitch error is equivalent to one 
degree of FD error. Thus, a six degree pitch limit would corres- 
pond to a limit of six degrees of FD error. Similarly, the 2.6 m/s 
airspeed limit would result in about a 5 degree error on the FD. 
Thus, if the pilot wishes to maintain pitch or airspeed within 
tolerable limits, he should keep the FD error below 5-6 degrees. 

At the same time control and control-rate should be maintained 
within the previously stated limits (Table 2) . If this logic is 
applied to both longitudinal director laws, the cost functional 
weightings are those given in Table 10. 

*And the fact that it works so v/ell in single-loop laboratory 
tasks . 
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Table 10 

LONGITUDINAL FLIGHT DIRECTOR WEIGHTINGS 


| Variable 

"Limit" 

Weighting 

FD 

s 

6 (deg) 

.028 

6 e 

9 (deg) 

.012 

6 e 

50 ( deg/s) 

.0004 

fd n 

1 ( "DOT" ) 

1 


29 (deg) 

.0012 

S N 

100 (deg/s) 

.0001 

- ■ • 


Analysis of the longitudinal flight director system was 
carried out using a cost functional corresponding to case (ii) 
("standard weightings") and a cost functional corresponding to 
case (iii) (Director-Weightings) . The results are compared below. 
The fourth strategy mentioned above was not investigated here, 
although it has substantial face-validi ty . The reason was that 
the corresponding cost-functional would be a combination (of sorts) 
of those for (ii) and (iii). Therefore, we would not expect sub- 
stantially different results from those we obtained. (Indeed, one 
might expect the results for (iv) to be intermediate between (ii) 
and (iii).) 
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Analysis and Results 

Comparison of Performance for Different Weightings. — Per- 
formance as a function of observation noise/signal ratio was computed 
for the 1% and 50% wind conditions and was averaqed across winds, 
as before. It was assumed that only flight director signals were 
displayed and that attention was divided eaually between stick and 
nozzle directors. The results for the two sets of weightings are 
presented in Figure 15. Also shown, is the result for the status 
display. It can be seen that in high workload situations (i.e., 
relative attention of 1/2 to 2) , performance with the status dis- 
play is nearly as good or better than performance with the flight 
director, depending on the choice of cost functional. In other 
words, if the pilot works hard enough, he can do about as well 
with the status display as with the flight directors. On the 
other hand, in the range of workload that is likely to be accep- 
table (1/2 to 1/8 relative attention) , the flight director system 
yields improved performance for a given workload (for both sets of 
weightings) . Moreover, the differences in overall performance 
obtained with the different weightings are relatively less signi- 
ficant . 


It is of interest to examine the nature of the differences 
between the results for the two sets of weightings. This is 
illustrated in Figure 16 where height and airspeed performance 
are compared for the two weighting conditions (for the 50%-wind) . 
It can be seen that the two cost-functionals lead to different 
performance tradeoffs. Height errors are regulated more tightly 
in the case where the director signals are weighted explicitly, 
but at the expense of relative degradation in airspeed regulation. 
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RELATIVE ATTENTION 


FIGURE 15 o Effect of Cost Function on Missed^Approach Probability 
with Interim Director (Averaged for all-winds) 
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FIGURE 16. Effect of Cost Function on Height and Airspeed 
Performance with Interim Director (50%-Wind) 
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For the remainder of this investigation we shall use the 
standard weightings that correspond to (ii) above. As noted 
above, we cannot decide a priori how the pilot will use the 
directors and, therefore, what the appropriate cost functional 
will be. The cost functional corresponding to the standard 
weightings appears to be a reasonable choice and allows us to 
make direct comparisons with no-director results. Moreover, the 
situation in which the pilot uses only director signals is some- 
what artificial; we suspect that when other instruments are 
available, the standard weightings are even more appropriate and, 
in addition, that the sensitivity to cost weightings is reduced. 

Allocation of Attention . When the fliqht director is 
presented along with the basic status display, it may be 
profitable to share attention among the various display elements; 
in any case, it seems likely that the pilot will do so. To 
explore this possibility, we conducted a model analysis to de- 
termine the optimal allocation of attention between the flight 
directors and the "standard" instruments of the status display. 

To simplify the analysis, we considered the directors and the 
status display as two entities between which total attention was 
to be shared. Then, the fraction of attention devoted to the 
flight director as a whole was allocated equally between the 
stick- and nozzle-director whereas the fraction devoted to the 
status display was split up among the height, pitch and airspeed 
in the same proportion as was found to be optimal for the no- 
director case (namely, .35, .1 and .5, respectively)* 
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Table 11 gives rms scores as a function of attention 
devoted to the flight director . Results are presented for a 
high-wind (1%), high-workload (P = -20dB) case; and for a median 
wind (50%), moderate workload {P Q = -14dB) condition. The scores 
indicate a good deal of insensitivity to allocation of attention, 
so that devoting anywhere from 25-75% attention to the flight 
directors is reasonable, with approximately 45% optimal. On the 
other hand, the results do confirm that it is undesirable to 
use only the flight directors. It is interesting to note that 
the penalty for paying full attention to either the flight 
director or the status is both poorer height reaulation and 
increased stick activity. Airspeed errors and nozzle-control 
activity are less sensitive to attention but do increase somewhat 
as attention to the flight directors increases. 
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Table 11 

EFFECTS OF ATTENTION SHARING BETWEEN LONGITUDINAL FLIGHT 
DIRECTOR AND STATUS DISPLAYS 





ATTENTION TO FLIGHT DIRECTOR 1 


VARIABLE 

0 

.25 

. 375 


. .5 

.75 

1.0" ' 

P =-20dB 
o 

a h (m) 

2.2 

H 

1.85 

1.82 

1.81 

1.84 

2.2 

1%-wind 

a (m/s) 
u 

1.33 

1.34 

1.34 

1.35 

1.35 

1.38 

1.46 


a Q (deg) 

1.51 

1.53 

1.55 

1.56 

1.57 

1.65 

1.94 


a Se (deg) 

1.85 

1.74 

1.77 

1.79 

1.81 

1.93 

2.42 


°6N (deg) 

12 

12.2 

12.2 

12.3 

12.3 

12.5 

13.0 


J 

1.31 

1.17 

1.16 

1.15 

| 1.16 

1.22 

1.57 

P =-14dB 
o 

a h (m) 

1.76 

1.25 

1.18 

1.15 

1.14 

1.14 

1.26 

50%-wind 

a u (m/s) 

.57 

.56 

.56 

.56 

.56 

. 58 

. 60 


a Q (deg) 

.78 

.76 

.71 

.72 

.72 

.77 

.87 


°fie (deg) 

.95 

.84 

i .86 

.87 

.88 

.97 

1.16 



4.33 

4.67 

4.71 

4.72 

4.73 

4.77 

4.78 


J 

.44 

. 30 

| .29 

.29 

.29 

. 30 

. 37 


4 












Report No. 2484 


Bolt Beranek and Newman Inc 


In Chapter 3, we suggested that the results for the idea- 
lized display configuration might be indicative of the maximum 
performance improvement to be expected from the flight directors. 
Moreover, a reasonable portion of that improvement would be due 
to display integration in the sense of a reduction in attention 
sharing. In Table 12 the scores for the idealized display are 
compared with those obtained (above) with the flight directors 
and with the flight director- status combination (assuming optimal 


Table 12 


COMPARISON OF RMS PERFORMANCE OF IDEALIZED DISPLAY 
AND FLIGHT DIRECTOR CONFIGURATIONS 
(1%-wind, P Q = -20dB) 


VARIABLE 

IDEALIZED 

FLIGHT DIRECTOR 

FLIGHT DIRECTOR + 

STATUS (OPTIMAL 
ALLOCATION OF ATTENTION) 

a h (m) 

1.53 

2.2 

1.82 

a u (m/s) 

1.30 

1.46 

1.35 

a Q (deg) 

1.77 

1.94 

1.56 

a 6e (dec,) 

2.0 

2.42 

1.79 

o {N (deg) 

11.9 

13. 

12.3 


allocation of attention) . These results demonstrate that, at 
least for the worst-case wind condition, the regulation of height 
errors when the flight director information is available is not as 
good as might be hoped for. Moreover, the fact that it is optimal 
to devote about half the attention to the status display, suggests 
that the director laws are not "integrating" the available infor- 
mation as effectively as possible. 
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Monitoring performance, as predicted by the ratio of rms 
estimation error to rms value, for various longitudinal display 
configurations is shown in Table 13. These results bear out and 
shed light on those of Tables 11 and 12. It can be seen that 
monitoring performance for none of the actual display configurations 

Table 13 


EFFECTS OF DISPLAY CONFIGURATION ON LONGITUDINAL 
MONITORING PERFORMANCE 
(1%-wind, P Q = -20dB) 



IDEALIZED 

urns 

FLIGHT DIRECTOR 

FLIGHT DIRECTOR 
PLUS STATUS 

1 °h /0 h 
1 e 

.11 

.42 

.30 

.36 

a h /a h 
1 e 

. 33 

.60 

.49 

.53 

°ee /o e 

.12 

.26 

.47 

.32 

°qe^°q 

.09 

. 16 

.24 

.19 

| 0 ue // °h 

L_i^ 

.59 

.78 

.62 


approaches that for the idealized display. The use of just the 
flight directors does lead to improved estimation of height errors 
but at the cost of a substantial penalty in pitch and airspeed 
estimation. When the status display is used in conjunction with the 
flight director a useful compromise among estimation errors for 
all variables is apparently achieved. Another interesting point 
is that estimation errors for height and sink-rate are less than 
for the status display. Apparently, the reduced thresholds and 
the zero-reference (lower residual noise) of the flight director 
offset the requirement for the processing of director signals to 
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obtain estimates of height and sink-rate. On the other hand, 
pitch, pitch-rate and airspeed are estimated better from the 
STOLAND status indicators (which are properly scaled and have 
references) than is possible via processing of the director 
signals . 

Performance and Workload . — It was just shown that 
improved performance may be obtained when the "interim" longi- 
tudinal director is added to the status display. Figure 17 
illustrates the average (over all winds) improvement in perfor- 
mance and the reduction in workload that may be expected for 
this display configuration. We assume that approximately 45% 
of the pilot's attention is devoted to the flight directors. 

The combination of status and flight director yields better 
window performance and a substantially reduced workload when 
compared with that attainable with either display alone. For 
example, if the criterion performance level is set at 1%, the 
workload index for the STATUS-FLIGHT DIRECTOR combination is 
approximately .45 whereas for either display alone it exceeds 
1 . 0 . 
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FIGURE 17. Effect of Longitudinal Displays on Missed Approach 
Probability 
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Lateral Director System 

The proposed "interim" lateral flight director system is 
described and discussed in [17], The approximation to that sys- 
tem used here is shown in Figure 18 (aloncr with the representa- 
tion of lateral dynamics) . It should be noted that we assume 
here that lateral flight path angle may be obtained directly 
rather than by means of the complementary filterincr technicrues 
of [17], Although this assumption is somewhat unrealistic, 
the idealization should provide a bound on the performance 
improvements that can be expected of the more practical system. 
The gains for the lateral director system are given in Table 14 
(they correspond to case 2F of [17] )» 

Vehicle dynamics for this analysis were the SAS-on, coordi- 
nated turn dynamics used previously (Appendix C) . Model para- 
meters were all set at nominal values. Thresholds and residual 
noises were assumed to be negligible for the lateral flight 
director. Cost weightings were the same as those used in the 
basic analysis (Table 2B) , so we assumed, effectively, the 
pilot's task is the same when the director is available as it 
is without it. 


Table 14 

GAINS FOR INTERIM LATERAL DIRECTOR 

.00638 volts/meter 
K c 1.55 volts/RAD 

1.61 volts/RAD 

D 

K 0 1.0 in displav/volt 

O 

T^ 1.0 sec. 
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Analysis and Results 

The lateral director system was analyzed in the same 
manner as was the longitudinal director svstem. Results are 
presented below. 

Allocation of Attention . — The effects on rms perfor- 
mance of sharing attention between the lateral flight director 

and the lateral status displays (considered as an entity) 
are shown in Table 15. It was assumed that the portion of atten- 
tion devoted to the status displays was allocated between the 
localizer and bank angle indicators in the approximately optimal 
3:1 ratio found in Chapter 3. 

Table 15 

EFFECTS ON LATERAL PERFORMANCE OF ATTENTION-SHARING 
BETWEEN FLIGHT DIRECTOR AND STATUS DISPLAY 
(1%-wind, P = -20dB) 


Attention to Flight Director 


Variable 

0 

.25 

. 5 

in 

• 

.8 

.9 

i- 1 

• 

o 

a y (m) 

5.1 

4.9 

4.8 

4.8 

4.7 

4.7 

4.6 

a* (m/s) 

2.0 

1.9 

1.9 

1.9 

1.9 

1.9 

1.9 









a ^ (deg ) 

5.0 

4.8 

4.7 

4.7 

4.7 

4.7 

4.7 

o, (dea) 

10. 6 

10.1 

9.9 

9.8 

9.8 

9.8 

9.9 

w 








J 

1.34 

1.22 

1.18 

1.16 

1.15 

1.15 

1.19 
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The results indicate that about 80-90% attention to the 
flight director is "optimal", but they are very insensitive to 
changes in attention. Referring to Table 8, column E (the 
idealized display condition) , we see that results for the 80% 
division of attention are cruite close to those for the idealized 
display; lateral error is about 3% greater for the flight director 
status combination and other variables are virtually identical. 
Even when only the flight director is available, there is not a 
significant increase in lateral error. Monitoring performance, 
not shown here, exhibited the expected trends (given in the 
earlier results). In general, then, the attention-sharing 
results indicate that the "interim" lateral flight director comes 
closer to achieving the improvements implicit in an idealized dis- 
play than does the longitudinal-director system. On the other 
hand, the improvements at this level of attention (-20dB) and wind 
condition are far from dramatic indicating that at high workload 
the status displays perform nearly as well as the director. 

Performance and Workload .— The improvement provided by 
adding the lateral flight director to the display at various 
"levels" of attention is shown in Figure 19. Averaged over all 
wind conditions, the probability of a missed approach without 
the flight director is 1.5 to 2 times greater than with it — 
at all levels of attention. Moreover, the improvement is 
greatest in the range of operation (attentions of . 5 to .125) 
that are likely to be most important. The situation for the 
50%-wind (also shown in Figure 19) is that the flight director 
provides even more substantial improvement. 

In terms of workload, an approach success probability of 
99% is unattainable, with reasonable workload, when all winds are 
considered; for the 50%-wind, the flight director reduces the 
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FIGURE 19. Effect of Lateral Director on Missed Approach 
Probability 
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Workload Index from about 1.4 to about .3. For the all-winds 
average, a 95% success probability requires a workload index 
of about .7 in the no-director case as opposed to about .25 when 
the director is available. In general, the curves of Figure 19 
indicate that addition of the flight director will reduce the 
lateral workload by a factor of 2-4 for success probabilities 
that can be achieved, again with greatest improvement in the 
range of most interest. 

Combined Longitudinal and Lateral Performance 

The total longitudinal-lateral approach task with the 
flight directors was analyzed in exactly the same fashion as was 
the EADI-status display configuration. First, the "optimal" 
allocation of attention between longitudinal and lateral control 
was determined, assuming attention sharing within tasks is the 
"optimal" for that task being performed alone. It was found 
that attention should be shared equally (approximately) between 
the longitudinal and lateral control tasks. Given this atten- 
tion-split the overall missed approach probabilities may be 
computed as functions of the relative attention devoted to the 
combined task. The results are presented in Figure 20. 

When the average of all-winds is considered, the addition 
of the flight-directors reduces the miss probability by about a 
factor of two, with the most improvement in the lower attention 
levels. Even greater improvement (4-7 times better) is evidenced 
for the 50%-wind condition. It may be seen in Figure 20a that the 
missed approaches, for the all-winds average, are due largely to 
the lateral task; although not shown, this is even more true for 
the 50%-wind condition. Table 16 illustrates the effect of the 
flight directors on workload. It can be seen that for the 50%-wind 


84 





1 


1 





Report No. 2484 


Bolt Beranek and Newman Inc. 


condition, the directors cut the workload by at least a third in 
the range of success probability of 95-99%. A similar reduction 
in workload is possible for the all-winds average, but the pro- 
bability of success is much reduced. (The workload index for a 
95% success probability is about .7 with the directors as opposed 
to 2 without them.) 


Table 16 

EFFECT OF INTERIM-DIRECTORS ON WORKLOAD 
( 50%-wind) 


Workload Index 


Probability of 
Missed Approach 

Without Director 

With Director 

1% 


.66 

2% 

1.6 

.47 

3% 

1.25 

.36 

4% 

1.04 

.3 

5% 

.88 

.28 


SUMMARY 

The addition of the "interim" longitudinal and lateral 
flight director systems improves performance considerably and 
reduces pilot workload by a significant amount. When the 
average of all-winds is considered, reducing the probability of 
a missed approach to 5% still requires a high workload. However, 
for the 50%-wind condition workload requirements seem well within 
capabilities . 
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Lateral offset errors are largely responsible for the 
missed approaches. However, the lateral director comes very close 
to matching the performance that might be obtained with an idea- 
lized display whereas the longitudinal director is less successful 
in this respect. In addition, the attention-sharinn studies 
suggest that the lateral director system is more effective in 
combining or "integrating" displav information than is the loncri- 
tudinal system. Thus, the poorer lateral performance would appear 
to be inherent in the dynamics, and further improvement in the 
lateral-directional axes via better display design mav be cruite 
difficult to achieve. 


87 



Report No. 2484 


Bolt Beranek and Nev/man Inc. 


88 



P^CEnTNTG PAGE BLANK NOT FILMED 


Report No. 2484 


Bolt Beranek and Newman Inc. 


AN APPROACH TO FLIGHT DIRECTOR DESIGN 

In the preceding chapter, proposed interim flight director 
laws were analyzed. These laws were derived on the basis of classical 
automatic and manual control theory in independent studies [17, 18]. 

We now describe and analyze a method for designing flight director 
laws that is based on, and is consistent with, the optimal control 
model that has been used throughout this investigation. 

The approach and the director laws explored here represent 
only a preliminary attempt at flight director design. Other ap- 
proaches based on the optimal control model are certainly possible. 
Moreover, a realistic flight director design procedure would involve 
consideration of factors (such as wind shears, beam capture) that 
were neglected in this analysis. Finally, we consider only longi- 
tudinal control. Nevertheless, we expect that the preliminary 
design developed here is sufficient to indicate the degree of 
performance improvement and workload reduction that is attainable 
in practice and to demonstrate the design approach. 

Director Design Philosophy 
Design Goals and Constraints 

The basic purpose of a flight director is to reduce the 
pilot's workload and/or to reduce system, errors. In terms of 
the optimal control model, this goal is accomplished by pro- 
viding the pilot with information necessary for flight control 
in such a way that the effects of display-related sources of 
pilot remnant are minimized. 


Preceding page blank 


89 



Report No. 2484 


Bolt Beranek and Newman Inc. 


The director variables must be composed of signals that 
can be generated by available aircraft measurement devices. These 
signals are commonly combined in such a way as to make the combined 
director-vehicle dynamics approximate a K/s-like behavior so that 
the need for pilot lead is minimized. Considerations of pilot 
acceptance also suggest that required control inputs and resultant 
vehicle motions be similar to those that are appropriate to flight 
with conventional displays. In addition, the director signal 
should provide the pilot with a good indication of instantaneous 
flight-path and attitude errors so that frequent reference to 
status displays is not reouired. 

If the director is a control director, a director signal 
must be generated for each control variable. It is likely that an 
ideal design should require little or no pilot coupling. That is, 
the output of a given director would command a control response 
along a single dimension. 

Preliminary Design Procedure 

The preliminary design procedure outlined below allows 
one to approach many of the above design goals in a relatively 
straightforward manner. The following design steps have been 
followed in deriving the control-director laws used in this study: 

1. Define the control situation in terms of system 
dynamics, input characteristics, sensory informa- 
tion, performance cost functional, and pilot 
parameters . 
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2. Use the pilot/vehicle model to predict the pilot- 
generated feedbacks between each display variable 
and each control variable. These feedbacks are 
designated as the "internal transfer functions" 
in the current implementation of the model. 


3. Approximate each of the internal transfers by a 
first- or second-order filter. (This approximation 
is done primarily to minimize the computational 
requirements for evaluating the proposed director 
design . ) 

4. The commanded control signal is generated by summing 
the outputs of the transfers between all sensor 
variables and the control variable appropriate to the 
director. The director signal is thus expressed as 


D i (s) = Zy . (s) • T' i . (s) 

-i J J 


where is the director signal appropriate to the i t ^ control 
variable, y^ is the j th sensor variable, and T'„ is the approxi- 
mate describing function between the sensor variable and the 

i th control variable. 


In essence, the flight director laws are designed to perform 
the equalization and cross-coupling that the pilot would other- 
wise have to do. With the director in the system then the 
pilot's task is basically that of generating a control response 
proportional to the deflection of the corresponding director 
indicator. Thus, cross-coupling should be at a minimum, and 
the pilot's response strategy should be approximately that of 
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a pure gain at low and mid frequencies. Since the transfers 
used in generating the director laws are only approximations to 
the predicted pilot transfers, the compensation provided by the 
flight director will be suboptimal with respect to the stated 
mission requirements. Consequently, a certain amount of cross- 
coupling is expected in the "pilot's" control strategy with the 
director. If the sensor variables and cost functional used in 
the design procedure are the same as would apply to the control 
task with a more conventional display panel, the characteristics 
of control and vehicle motions should not be appreciably changed 
by the use of the director. Improvement in performance and work- 
load reduction will accrue from the reduced effects of pilot- 
induced noise. If additional sensory information (such as linear 
and/or rotational accelerations) is used in generating the director 
signal, further improvement may be expected. 

Definition of the Control Situation 

The first step of the design procedure calls for a definition 
of the control situation. A description of the vehicle dynamics, 
input characteristics, and to a large extent sensory information, 
are known quantities and can be readily specified in a format 
appropriate to the pilot/vehicle model. The selection of approp- 
riate cost functional and, in this case, of the model noise/signal 
ratios reauires some careful thought. The choice of cost functional 
for analysis of flight director performance was discussed in 
Chapter 4. For the reasons given there (and with the same caveats), 
we use in this analysis the same cost functional that was employed 
in the no-director studies. 
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The selection of model noise/signal ratios for flight 
director design depends on whether one views the benefit of the 
director as primarily the reduction of system errors or the 
reduction of pilot workload. If the pilot is expected to maintain 
a high level of workload so that he can minimize errors, the noise/ 
signal ratios used in the analysis should be those appropriate to 
maximal effort. (Observation noise/signal ratios of -20dB and 
motor noise/signal ratios of -25dB are nominal values.) On the 
other hand, if the director is intended mainly to allow the pilot 
to maintain performance with reduced workload, then the director 
should be optimized for substantially larger noise/signal ratios. 

Although a director designed for a low-noise situation will 
allow better performance under conditions of maximal effort, the 
alternative design will be less sensitive to pilot noise and 
should thus be more "forgiving" of non-optimal pilot behavior. 
Accordingly, the director laws designed and evaluated in this 
section have been obtained by computing predicted pilot describing 
functions for a high-noise situation. 


Derivation of Director Laws 

Pilot describing functions were obtained from a steady- 
state analysis of pilot/vehicle performance under conditions 
of high noise/signal ratios. Vehicle dynamics and input charac- 
teristics were the same as those used in previous analysis and 
the same cost functional was used to describe mission requirements 
(Table 2A) . Wind gust intensities corresponded to the "1-percent" 
wind condition. The sensory variables used in designing the director 
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laws were the same as those assumed to be obtained from the STOLAND 
display (glide path error and rate, pitch error and rate, and 
airspeed error) . No additional acceleration inputs were considered. 


A time delay of 0.2 seconds was assumed for the pilot, and 
both observation and motor noise/signal ratios were set at -lOdB . 
Attention-sharing was assumed not to be a factor, and thresholds 
and residual noise sources were ignored. 


The internal pilot describing functions (magnitudes only) 
are shown for the elevator and nozzle controls in Figures 21 and 
22, respectively. In the process of deriving the flight-director 
laws, these curves were approximated by transfer functions of 
second-order, critically-damped, low-pass filters. In order to 
minimize the number of state variables needed to describe the 
director characteristics, the critical frequencies of all 
responses corresponding to a given control variable were made 
identical. Thus, each of the two director signals was represented 
as follows: 


D i (s) 




^1+s/wjJ 


2 

y j (s) 


where uk is the critical frequency of the filter associated with 
the i th control variable and ^ is asymptotic low-frequency 
behavior of the approximate transfer function relating the i ^ 
control variable to the j ^ sensory input. 


Second-order approximations to the predicted internal 
describing functions were obtained by visual inspection. The 
resulting director parameters are shown in Table 17. The units 
of the low-frequency gains are in terms of relevant display and 
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21. Magnitudes of Internal Describing Functions for Elevator Control 
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control variables. For example, the gain associated with the 
contribution of height information to the elevator director has 
units of degrees (of control surface deflection) per meter (of 
height error) . Critical frequencies are in radians/second. In 
the subsequent discussion, we shall refer to the directors using 
these parameters and laws as "model-based directors". 

Table 17 

PARAMETERS FOR MODEL-BASED LONGITUDINAL DIRECTOR LAWS 


Director 

Sensory 

Variable 

Critical 

Frequency 

Low-Frequency 

Gain 


Height 


0.25 


Sink Rate 


1.8 

Elevator 

Pitch 

5.0 

0.9 


Pitch Rate 


0.8 


Airspeed 


0.5 


Height 


1.1 


Sink Rate 


6.3 

Nozzle 

Pitch 

3.5 

0.8 


Pitch Rate 


1.0 


Airspeed 


6.3 
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Analysi s and Resul ts 

As opposed to the design phase, motor noise/signal ratios 
were fixed at approximately -25dB, and observation noise/signal 
ratio was a parameter for the analysis. Attention was assumed 
to be shared equally between the elevator and nozzle director 
signals. (No attempt was made to find the optimal allocation 
of attention. On the basis of previous analysis, it was 
assumed that performance would be relatively insensitive to 
attentional allocation.) Thresholds for the stick and nozzle 
directors were set at the same values that were used for the 
analysis of the interim directors. Other model parameters 
were the same as those used in the design of the laws for the 
directors . 

Performance . — Predicted rms performance scores for the 

1%-wind condition and an observation noise/signal ratio of -20dB 
are shown in Table 18. Also shown for comparison are the scores 
predicted for the status display without a flight director 
(Condition A of Figure 8) the scores associated with the idea- 
lized display condition (Condition F of Figure 8) , and the 
scores for the interim director. The interim director 
results are presented as a matter of interest and comparison 
of them with model-based director results is not wholly 
warranted. The interim director design and analysis were 
based on slightly different dynamics (See Appendix C) . Further- 
more, the interim director design did consider a number of 
factors other than gust regulation. On the other hand, compa- 
rison of model-based director results with those of the basic 
status- and idealized -display conditions is entirely appropriate. 
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The performance variable most effected by display para- 
meters is the rms height error, as expected from previous results. 
The score predicted with the flight director is about 26% less 
than the score predicted for the status display. (The idealized 
display yields about a 35% reduction with respect to the status 
display.) A similar reduction is predicted for the sink rate 
error score. A reduction of about 10% is predicted for rms 
stick and stick rate. Other performance scores are virtually 
unchanged . 


Table 18 

COMPARISON OF RMS PERFORMANCE WITH 
VARIOUS DISPLAY CONFIGURATIONS 

(1%-wind, P Q = -20dB) 
Display Condition 


Variable 

Without FI 
Status 

ght Director 
Idealized 

With FI. 
Interim 

Lght Director 
Model-Based 

o h (m) 

2.3 

1.5 

2.2 

1.7 

(m/s) 

.68 

.55 

.72 

.55 

a Q (deg) 

1.8 

1.8 

1.9 

1.8 

(deg/s) 

1.2 

1.2 

1.6 

1.1 

a u (m/s) 

1.3 

1.3 

1.46 

1.3 

a. (deg) 
°e 

2.1 

2.0 

2.4 

1.9 

at, (deg/s) 
°e 

7.4 

7.0 

8.9 

6.6 

a {N <deg) 

12. 

12. 

13. 

12. 

at (deg/s) 
°N 

21. 

9 

i 

21. 

23. 

21. 
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Except for improved flight-path performance, then, vehicle 
motions and control responses are essentially the same with 
and without the model-based flight director. 

Estimation or monitoring performance for the above case 
is quite interesting and is shown in Table 19. It may be seen 
that height estimation-errors are greater percentage-wise than 
for any of the other configurations. Pitch, pitch-rate and 
airspeed estimation performance for the model-based director 
is about the same as for the status display and better than for the 
interim director; sink-rate errors are smaller than for 
obtained with STOLAND and the same as for the interim director. 
Thus, when compared with the status display, the improved height- 
regulation evidenced in Table 18 appears to be the result of 
better sink-rate estimation (even though height errors them- 
selves are not so well estimated). On the other hand, the 
improvement in performance over the interim director seems 
to result from better estimation of pitch and airspeed. 

Table 19 

MONITORING PERFORMANCE FOR MODEL-BASED 


DIRECTOR SYSTEM 
(1%-wind , P q = -20dB) 

Display 


Variable 

Idealized 

Status 

"Interim " 
Director 

"Model-Based" 

Director 

°h e /0 h 

.11 

.42 

.30 

.56 

C h e /0 h 

.33 

. 60 

.49 

.48 

‘V’e 

.12 

.26 

.47 

.28 

°q. /0 <S 

.09 

.16 

.24 

.18 

a u e /o u 

.54 

.59 

• 

00 

.59 
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Allocation of Attention . — To assess the extent to which 
the flight director is providing the required "mix" of system 
variables, the effects of sharing attention between the model- 
based director and the status display were investigated (Table 20) .* 
The "optimal" allocation to the flight director is to devote 
80-90% attention (recall, for the interim directors about 45% 
to the directors was optimal). Again, the results are highly 
insensitive to the exact allocation. However, unlike the case 
of the interim director, in this case devoting full attention 
to the flight directors imposes virtually no penalty. Thus, 
it appears that the model-based director laws provide signals 
that are an effective mix of the system variables. 

Table 20 

EFFECTS OF ATTENTION-SHARING ON PERFORMANCE 
WITH MODEL-BASED. FLIGHT DIRECTOR 

(1%-wind, P q = -20dB) 



Attention to Fliaht Directors j 

Variable 

0 

.4 

. 5 

.75 

. 8 

.9 

1.0 









a h (m) 

2.2 

1.77 

1.75 

1.72 

1.72 

1.72 

1.75 

a u (m/s) 

1.33 

1.32 

1 . 32 

1.32 

1 . 32 

1.32 

1.32 

a e (deg) 

1.51 

1.76 

1.76 

1.76 

1.76 

1.76 

1.76 

a 6 e (deg ) 

1.85 

1.99 

1.98 

1.96 

1.96 

1.95 

1.94 

Ot (deg) 
6 N 

12 . 0 

12.0 

12.0 

12.1 

12.1 

12.2 

12.3 

J 

1.31 

1.14 

1.13 



1.12 

1.11 

1.10 



1.12 


*As in the earlier studies of Section 4 , it was assumed that 
the flight directors and status displays could be considered as 
entities, with attention within these entities divided as before. 
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Work! oad . — In order to assess the degree to which work- 

load can be reduced by the use of the flight director, we examine 
the relationship between the probability of a "missed approach" 
and "attention". A missed approach is defined, as before, as 
the situation where either height and airspeed errors exceed 
their respective "limits" of 3.7 meters and 2.6 meters/second. 
Attention is related inversely to the observation noise/signal 
ratio, with a relative attention of unity associated with a 
ratio of -20dB. 

The relation between predicted performance and attention 
is shown for the 1%-wind condition in Figure 23a. In addition 
to results for the model-based director, curves for the status 
display and the interim director plus status display combi- 
nation* are shown for comparison. For the 1%-wind condition 
performance is still poor for the model-based directors, but 
it is appreciably better than for either of the other display 
configurations. In particular, the model-based flight director 
reduces significantly the sensitivity of performance to observa- 
tion noise (both display-related and human related) and, there- 
fore, shows relatively greater improvement at lower levels of 
pilot attention. This is very evident in the 50%-wind case (not 
shown) where, with the model-based director, the probability 
of a missed approach is essentially zero throughout the range 
of attention-levels investigated. This insensitivity was, of 
course, a prime objective of the design approach. 


* Recall that performance and workload for the combination was 
better than that for the interim director alone. 
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The effect of attention on performance averaged across 
gust conditions is shown in Figure 23b. The value of the 
flight director is now even more apparent. The predicted 
workload requirements are substantially lessened even for 
relatively stiff performance demands. For example, if we 
require a 99-percent probability of a successful approach, a 
relative pilot attention of slightly greater than unity is 
required when no director is provided. The interim director- 
status combination reduces this " attentional demand” by about 
a factor of 2 (i.e., to .4 - .5). With the model-based direc- 
tor, however, attention requirements are reduced by about a 
factor of 10. Conversely, performance is improved for a pilot 
operating at a constant level of attention. For a relative 
attention of 1/4 the predicted probability of a missed approach 
is reduced from about 5 percent for the status display and 2 
percent for the interim-director-status combination to 
around .6 percent for the model-based director. 

Predicted Pilot Describing Functions . — One of the 
design goals set forth earlier was that the flight director 
should allow the pilot to adopt a control strategy that re- 
sembles a simple gain at low and mid frequencies. We pre- 
dicted that the design procedure adopted in this study would 
meet this reauirerent by allowing the director laws to perform 
the required equalization. It was also anticipated that cross- 
coupling in the pilot's response strategy would be unnecessary 
with a properly designed set of flight directors. 

Inasmuch as the model for the pilot is relatively "free 
form", pilot transfers will in general be predicted between 
all display and all control variables. Thus, for the control 
situation investiaated here, there are two sets of predicted 
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FIGURE 23. Effect of Attention on Missed Longitudinal Approach 
for "Model-Based" Directors 
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pilot describing functions to consider: the "direct" transfers 

which relate each control response to the corresponding director 
command, and the "cross" transfers which relate control responses 
to commands on non-associated directors. In cases where pilot 
cross-coupling is unimportant, the magnitudes of the predicted 
cross transfers should be numerically small. 

The predicted direct transfers are shown in Figure 24. As 
expected, these transfers approximate a pure gain at freguencies 
up to about 4 rad/sec (which is beyond gain-crossover for flight- 
path and attitude control) . The high-f reguency peaks in the 
amplitude ratios are typical of actual pilot response behavior 
obtained in K/s tracking situations. 

Predicted cross transfers are shown in Figure 25. The 
freguency-dependency of the phase-shift indicates that both 
describing functions are non-minimum-phase (i.e. , there 
are zeros which have positive real parts) . 

In order to determine whether or not the magnitude of the 
cross transfers are small enough to be neglected, we must compare 
the open-loop describing functions for the direct and cross 
paths. (The open-loop transfer is defined as the cascade combi- 
nation of a predicted pilot describing function and the associa- 
ted vehicle transfer function.) Such a comparison (not shown 
graphically in this report) reveals that the magnitudes of the 
cross control paths are not substantially less than the magni- 
tudes of the direct paths at all freguencies. Thus, we cannot 
claim that the predicted pilot cross couplings are numerically 
small . 
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FREQUENCY (RADIANS/SECONDS) 


FIGURE 24. Predicted Director-Control Describing Functions: 
"Direct" Transfers 
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FREQUENCY (RADIANS/SECONDS) 

FIGURE 25. Predicted Director-Control Describina Functions: 
"Cross" Transfers 
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A true test of the importance of cross-coupling would be 
to determine the levels of performance and workload that would 
be obtained if cross coupling were prohibited. There is no 
simple way to make this test at present, however, because 
current implementation of the pilot/vehicle model does not 
allow for such a constraint on the predicted control strategy. 
Thus, it is presently not possible to determine conclusively 
whether or not the pilot must introduce cross-coupling in order 
to realize the benefits of the flight director. 

SUMMARY 


In this section, an approach to designing flight director 
laws based on the "optimal-control model" of the human operator 
was suggested. Director laws for longitudinal control were 
developed using this approach. These were evaluated in the 
same fashion as was the status display and interim-director 
analyzed previously. It was found that the model-based 
director provided improved performance at substantially re- 
duced workloads and thus achieved its major design objectives. 
On the other hand, the results do not substantiate the belief 
that the need for control cross-coupling would be reduced by 
this design procedure. Further work is necessary to evaluate 
this aspect of the design. 

It should be reemphasized that the design procedure 
presented here is only in a preliminary stage. Consideration 
was not given to aspects of the design, other than gust regula- 
tion, and design compromises that are perhaps inevitable in 
practice were unnecessary here. Nonetheless, we believe the 
results to be highly encouraging and worthy of further inves- 
tigations . 
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RESPONSE TO WIND SHEARS 

In previous chapters, we have examined the performance of 
proposed vertical situation displavs in crust-regulation tasks. 
Steady-state analysis was used to evaluate performance and work- 
load at the approach window. However, the pilot's control task 
involves more than gust-recrulation. An important aspect of his 
task is compensation for errors introduced by winds with a non- 
zero mean component. These "mean-winds" will, in general, vary 
with altitude. The rate of variation of mean-wind speed with 
altitude is referred to as the shear variation. In the sequel, 
we shall often refer to these altitude-dependent winds as wind- 
shears . 


The mean-wind may, itself, be described in statistical 
terms, i.e., the wind direction and speed is a random variable. 
However, in a given approach- to-land ing , a specific "sample" 
mean-wind is encountered. It is the response to some particu- 
lar samples that we shall be concerned with here, rather than 
to the distribution as a whole. 

As noted in Chapter 2, the modification of the optimal 
control model for the human operator to account for pilot 
adaptation to disturbances with a " time- varying" mean is 
described in Appendix A and the equations needed to compute 
the response to a particular "sample" disturbance are given 
in Appendix B. It is noteworthv, that the computed model- 
response to a sample-disturbance is a random variable because 
of the inherent randomness of the human. In other words, the 
model predicts the distribution of responses that would be 
obtained if the approach were repeated by the pilot a number 
of times (theoretically, an infinite number) with the same 
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"sample" disturbance. For this analysis, we shall assume that 
the system disturbances will consist of "sample"-shears and 
turbulence that is described as in Eouation (12) , Chapter 3. 
Therefore, there will also be a random component of the response 
that is attributable to the statistical model for turbulence. 

Longitudinal Analysis 

Modelling the Wind-Shears . — In modelling the wind-shear 
effects there are both dynamic and kinematic effects to consider. 
In addition, from the standpoint of implementation of the optimal 
control model, it is desirable (though not necessary) to convert 
the altitude dependence of the shear to an " eauivalent" time 
dependence. In this section we discuss the approach taken to 
these issues concerning modelling of the wind-shears. We consi- 
der only the horizontal wind-shears. 

Figure 26 illustrates the pertinent geometry. The air- 
craft's altitude (h) is given bv 

h = h + <5h = R tan r + 6h (16) 

n o 

where h is the "nominal" altitude, i.e., the altitude of the 
n 

glide-slope at the aircraft's range, P, and Sh is the altitude 
error. The rate of change of the nominal altitude may be 
expressed in terms of the ground speed (or range-rate) . 

h = R tan T z (u + u) tan T (17) 

n o o o 
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where 

u = perturbation in around speed 

U q = nominal airspeed 

The aircraft's sink-rate is 

h = V sin y ~ (L T o + u) sin ( r Q + Ay) 

s (U + u) (sin r + cos r • Ay) 
o o o 

and 

6h = h - h 

n 

X (U + u) cos r o 0 Ay = (U Q + u) cos r o • (0-a) 

or 

cos r o ° w (18) 

Let u and w be the zero-wind lonqitudinal and vertical airspeed 
a a 

perturbation components and 

) (19) 


u = u + u Mt . 
a MW 


w = w + W MI1 
a Mw 


6h = (U + u) cos f o " 0 
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Then, the 

dynamic 

effects of 

the 

shears are 

accounted 

for by 

modifying 

the 

eauations 

of 

motion 

(A.ppendix 

C) 

, ★ 

to 

\ 


• 

u = 

X 

u 

(u - 

U MW 

+ 

V 

+ 

X w 

(w - 

W MW + 

w ) 
g 

+ • • . 


( 1 -Z • ) w = 
w 

Z 

u 

(u - 

U MW 

+ 

V 

+ 

Z 

w 

(w - 

W MW + 

w ) 
a 

+ ... 

1 

; 

> ( 20 ) 

• 

a = 

M 

u 

(u - 

U MW 

+ 

V 

+ 

M 

V! 

(w - 

W MW + 

V 

+ ... 

i 

> 



Eauation (18) is used to account for kinematic effects (in the 
sense that the eauation for altitude errors in Appendix C is no 
longer valid) . However, Eauation (18) is nonlinear if u is sub- 
stantial, as is the case for the winds to be considered here. 

Rather than neglect the nonlinear terms in (18) we shall use the 
"average" wind velocity during approach in computing glide-path 
errors. Thus, if u.„ , (0) is the initial wind velocity and au,_ t 
is the change in wind velocity during the approach, we let 

_ \ 

5h = U cos F • 6 - £ 7 — cos r • w 
o l< o 

o 

U = u 0 + U MU (0) + Au MW /2 

This choice tends to minimize the maximum error associated with 
assuming constant ground speed for linearization purposes. 

The wind-shears to be considered here are enumerated in 
Table 21. These winds are idealizations of more exact models 
for mean-winds. They were used in this analvsis so as to be 
compatible with a concurrent simulation study at Ames Research 
Center. We now show how these winds may be represented as time- 
varying disturbances; it turns out that this can be done with 
considerable fidelity. 

*x . X Z ! Z Z • , M and M are dimensional stability derivatives, 
u' w' u' w' w' u w 
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Let be the longitudinal fralong track) wind component of 
interest (crosswinds may be treated analogously) . Given the 
profiles of Table 21 , we may write 


U 


MW 


= U 


MWO 


+ a h 


( 22 ) 


where a is the change in windspeed with altitude, i.e., the 
shear- variation. Thus, using (17) 


= a h a a h = a (U + u) tan r (23) 

MW no o 

Dif ferentrating (23) gives 

U mT = (a tan V ) u ( 24 ) 

MW o 

This equation along with the dynamical equation for u (Eouation 
(20)) allows the wind shear to be expressed in terms of other, 
non-input related, state-variables. An even simpler represen- 
tation is possible and seems warranted in terms of the objecti- 
ves of this study. If we assume that the pilot attempts to 
maintain airspeed (u s0) , then Equation (19) implies that we may 
substitute u MW for u in Equations (23) and (24) . But 


u. __ ~ U._, cos 0 
MW MW o 


(25) 


and, the state-variable representation for the mean wind is 


X 1 U MW x 2 ' x l (0) U MW (0) 


x„ - = (a tan T cos 0„)U„„ 

2 MW o o MW 


> (26) 


= (a tan r Q cos e Q )x 2 ; x 2 = a(U Q + u MWQ )tan r o 

so 

where x 2 is the mean-wind velocity at the onset of the shear, 
so 
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In Equation (20), we set w^ = sin 0 , completinq the 

model for horizontal wind-shears. 

Table 21 

WIND-SHEARS FOR LONGITUDINAL AND LATERAL ANALYSIS 


Wind 

BESSsB ' 1 

Initial 

Speed 

Final 

Speed 

Final 

Altitude 

Decreasing 

Tailwind 

152m (500 ft) 

15.45^(30KTS) 

5.15^(10KTS) 

0 

Increasing 

Tailwind 

152m 

-5.15^(10KTS) 

s 

+5.15- 

s 

0 

Increasino 

Crosswind 

0 

152m 

-2. 575|(5KTS) 

+7.725^(15KTS) 

0 


+ Indicates tailwind or crosswind from left side 
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Scenario . — An analysis of the wind-shear performance 
as extensive as that conducted for the steady-state conditions 
was not performed. Instead, some typical results for the 
status display configuration were obtained. A major purpose 
was to see how performance with wind shears compared to that 
for the steady-state analysis. A second important objective 
was to test the sensitivity of the model results to various 
alternative assumptions . 


The basic scenario for the longitudinal analysis involved 
starting at an initial range of 1500m with a constant wind velo- 
city corresponding to the value at h = 152m (R c 1160m) . In 
all cases, turbulence having the spectral form of Equation 12 
was employed. The scale-lenqths were not varied with altitude; 
they were set at the constant value appropriate to the decision 

height. Gust intensities were a = 1 . 5m/s , = ,58m/s, the 

u q d 

10% wind condition. In all cases where the STOLAND displav was 
used the glide-path error shifted from an angular presentation 
to a height presentation at R s 575m (h ~ 75m t 250 ft.) , as 
indicated in [7] . Thresholds and residual noises were adjusted 
accordingly. Attention among display variables was assumed to 
be allocated in the manner that was found to be optimal for steady- 
state (Table 3) . The above-mentioned conditions were fixed 
throughout the analysis. Several variations with respect to 
other factors were investigated. These variations are described 
below. 


The aircraft was assumed to be on the glide slope, but two 
types of initial trim were investigated. 
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Trim Condition A . — The component of the wind along the 
aircrafts x-body axis is trimmed out but the z-bodv axis component 
is not, so w(0) = 0. From Ecuation (21) it can be seen that this 
implies 6h(0) = 0. However, there will be, for this condition, 
an initial acceleration along the z-body axis. 

Trim Condition R . — Both x and z components of the ini- 
tial wind are trimmed out. This implies an initial w(0) and a 
corresponding initial glide path error-rate. For the decreasing 

tailwind, this rate is 6h{0) = 1.28m/s whereas for the increa- 

0 

sing tailwind, <5h{0) =-.31 m/s. 

It had been noted in [ 6 ] that the pilot might vary his 
gains so as to tighten control as the decision-height was 
approached. An analogous result can be obtained with the model 
by making cost functional weightings range-dependent. Indeed, 
if it is assumed that the pilot attaches a fixed penaltv to 
angular deviations from the glide-path, rather than linear 
deviations, then the weighting on height-errors will be range- 
dependent. We investigated two conditions with respect to this 
"aain-scheduling " . 

Constant Gains : The cost functional was fixed and eoual 

to that used in the steady-state analysis. 

Varying Cains ; It was assumed that the annular glide- 
path error corresponding to the "window" dimension was 
appropriate. Then, the qains were chanced in three 
stages. This was accomplished by changing the weighting 
on linear-height errors according to the following 
schedule: 
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1500m 

< 

R 

< 

1160m 

• 

• 

Q h = .0029, 

q^ = .0326 

116 0m 

< 

R 

< 

575m 

: 

o h = .0117, 

q^ = .133 

57 5m 

< 

R 

< 

230m 

; 

q^ = .073, 

% m - 83 


Thus, the weighting over a range-interval corresponded to the 
weighting appropriate to the end-point of that interval, a con- 
servative choice. The intervals were chosen, as a matter of 
convenience, so that the end-points corresponded to points 
where other changes were required. 

It was expected that the human 1 s time delay would increase 
scores hut would not alter the basic character of the results. 

Because inclusion of the time delay increases significantly the 
costs of the time-varying computation, we decided to assume the 
time delay was zero. However, a comparison case in which a time 
delay of .2 sec. was used was run to illustrate the differences 
one might expect from includinq time delay. 

R e s u 1 1 s . — Mean and standard deviation scores, at the 

decision-height, are compared in Table 22. Several points are 
worth mention. First, the constant gain-no shear, zero-delay results 
are virtually equal to those of the corresponding steady-state 
analysis. This is more than a check on the program; it shows 
that in the absence of shears, the approximately 1250-1300m 
approach distance is sufficient for the errors to reach steady- 
state. Second, the effect of the wind-shear is more than just 
a non-zero mean response. It may be seen that the standard 
deviation of the tracking errors and of the controls is increased. 
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Table 22 

PERFORMANCE AT DECISION-HEIGHT FOR 
VARIOUS ANALYSIS CONDITIONS 

(Status Display, 0 = 1.5m/s, P = -20dR) 

o 



Steady-State 

Constant-Gains 1 

Vary ina -Gains 

Variable 

T = . 0 

x = .2 

T = 0 

T = 0 

T = .2 

T = 0 




No- 

Shear 

With 

Shear 



h (m) 

0 

0 

0 

• 1 6 

.22 

.31 

°h (m) 

1.73 

1.83 

1.72 

2.02 

2.2 

2.03 

h (m/s) 

0 

0 

0 

-.01 

-.016 

-.05 

(m/s) 

.48 

.50 

.47 

.63 

.70 

.63 

6(m) 

0 

0 

0 

-.6 

- . 6 

-.62 

o 0 (m/s) 

1.24 

1.29 

1.23 

1.27 

1.35 

1.26 

' 

u (m/s) 

0 

0 

0 

-.10 

-.095 

-.10 

o u (m/s) 

.94 

.96 

.96 

.98 

1.03 

1.00 

7 (deg) 

0 

0 

0 

-1.44 

-1.43 

-1.46 

o. (deg) 
c 

1.38 

1. 5 

1 . 39 

1.75 

1.94 

1.82 

Z (deg) 

0 

0 

0 

16.3 

16.2 

16 . 5 

c. (deg) 

N 

8.6 

8.5 

8.6 

11.5 

11.8 

11.7 
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This is a result of the coupling in the model of mean-and 
variance-responses that arises from the dependence of the obser- 
vation noises on the rms signal values and that of the motor- 
noise on rms control. In terms of missed approach probabilities, 
the increase in variance is the more significant effect. 

Third, the effect of time delay is, as expected, to increase mean 
and standard deviation of the error. The magnitude of the effect 
is largest for height-error with approximately a 35% increase in 
mean and a 10% increase in standard deviation. 

The final effect illustrated in Table 22 is that resulting 
from allowing the gains to varv. When compared with the constant 
gain case, it is seen that the principal effect at the window is 
on the mean-response. This effect, though large percentage-wise, 
is virtually negligible in terms of the missed-approach probability. 
The differences between constant and varying gains are more pro- 
nounced in the time-histories shown in Figure 27. These time- 
histories are curves passed through data points obtained every 
50m. The jump-discontinuities for the varying-gain case arise 
from the instantaneous gain-change and the associated jump in 
control value. These "jumps" apparently decay very rapidly. 

Because height errors are weighted less, thev are allowed to 
build up to a greater extent in the varying-gain case; however, 
as the threshold is approached the errors begin to be reduced 
rapidly (because of the higher weighting) , so that window per- 
formance is not significantly different for the varying- and 
constant-gain cases (Table 22) . Apart from differences in height 
control, the principal difference between the two-cases is in the 
initial transient in elevator and pitch. It seems clear that the 
early reduction in height errors for the constant-gain case is a 
result of a rapid pitch-down. 
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(29S) { 185) (IS) (29S) (18S) (7S) 

RANGE (m) 

(TIME-TO-DECtSiON HEIGHT, S) 


FIGURE 27. Effect of Varying Gains on Approach Trajectories 
a) Height anc! Nozzle Response 
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CONSTANT GAINS VARYING GAINS 



(29S) (18S) (7S) (29S) (18S) (7S) 

RANGE(m) 

(TiME-TO-DECISION HEIGHT, S) 

FIGURE 27. Effect of Varying Gains on Approach Trajectories 
b) Pitch, Airspeed and Elevator Response 


122 



Report No. 2484 


Bolt Beranek and Newman Inc. 


The above results were all obtained for trim-condition A (zero 
initial glide-path error-rate) . The only significant effect of 
starting in trim-condition B was on mean height error, namelv the 
mean-height error starting in condition P was about twice as 
large as that for starting in A. throughout the approach. The 
problem is, simply, that wi th the more adverse initial condition, 
there is insufficient time (on the average) to compensate for the 
constant wind before the wind velocity starts to change. Thus, 
at the onset of the shear (R s 1160m) the height error for B is 
about twice that for A, and it remains so. In terms of "window" 
performance, the effect is to increase the probabilitv of exceeding 
the window by about 1% . 

The excellent "window" performance obtained in the above 
analyses is somewhat misleading. As can be seen from Figure 26, 
the nozzle limit of 29° is less than one standard deviation from 
the mean for much of the approach (after the window-veloci ty 
starts changing). Thus, in a high percentage of the time the 
nozzle exceeds its limit. What this means is that the rate of 
descent capability of the aircraft, with throttle fixed , is in- 
sufficient for this wind. Further, the wind is of sufficient 
severitv to place the entire linearized analvsis in Question. 

On the other hand, the analvsis suanests that suitablv scaled- 
down winds may be adeauatelv controlled by nozzle and elevator 
inputs alone. 

In an attempt to get some estimate of the control- limited 
performance for the decreasing tailwind, a trajectorv was obtained 
for a case in which nozzle control and control-rate were heavily 
penalized in the reaion where excessive nozzle-control had been 
observed, i.e., in the 1160m < R < 575m interval. (Weightings 
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on nozzle and nozzle-rate were multiplied by 50) . To allow 
transient effects resulting from the initial constant wind to 
die out, the approach was started at 2000m. The result for 
height-error and nozzle-position is shown in Figure 28. It 
can be seen that nozzle responses to the shear variation in the 
heavily penalized region are virtually nil and the height errors 
increase accordingly. When the penalty is reduced, R < 575m, a 
relatively large mean-nozzle motion ensues in an attempt to reduce 
the mean-error. While some reduction occurs, the mean height error 
at the decision height is still three times the allowable error. 
Although these results are not intended to be definitive, they do 
illustrate the problem posed by this wind, when throttle is fixed. 

A constant-gain trajectory for the increasing tailwind was 

* 

also obtained and the results are shown in Figure 29. For this 
case both longitudinal and normal components of the wind were 
initially trimmed out (analogous to condition B above) . The 
window performance for this wind is compared with that for the 
decreasing tailwind in Table 23. Note that the turbulence 
intensity and spectrum is the same for the two cases. It may 
be seen that height errors are controlled more effectively for 
the increasing tailwind; airspeed is less-well controlled. The 
overall effect is a definite improvement. Two other points are 
worth noting. Referring to Table 23, we find that the standard 
deviation of the height and sink-rate errors for the increasing 
tailwind are very close to those obtained for the steady-state 
(no-delay) case. Thus, it appears that with the tailwind (which 


*As can be seen the wind approximation is not as close to the 
idealized wind as for the previous case, but certainly good 
enough. 
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FIGURE 
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8. Nozzle-Limited Response For Decreasing Tailwind 


125 


S821 



RANGE (m) 

(TIME-TO-DECISION HEIGHT, S) 


FIGURE 29. Response History for Increasing Tailwind Shear 
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Table 23 


COMPARISON OF WINDOW PERFORMANCE FOR 
DIFFERENT TAILWINDS 


Variable 

Decreasing Tailwind 

Increasinq Tailwind 

h(m) 

.16 

-.12 

o h (m) 

2.02 

1.69 

h (m/s ) 

-.01 

.006 

(m/s) 

.63 

.46 

0 (deq ) 

-.6 

-.49 

a Q (dec) 

1.27 

1.13 

u (m/s) 

-.10 

-.19 

• 

a u (m/s) 

.98 

.99 

IgCdeg) 

-1.44 

. 56 

o 5 (deq) 
e 

1.75 

1.3 

(deq) 

16.3 

-7.4 

a. (deq) 
°N 

11.5 

10.0 

Time-for- 
Approach ( s) 

~ 29. 

~ 43 
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starts out as a headwind) , enough time is added to the approach 
to allow the "pilot" (model) to reduce the errors to values 
commensurate with an approach of infinite length. The second 
point is that the nozzle-control requirements are not so excessive 
(in relation to capability) as for the decreasing tailwind. Thus, 
one might expect these results to correspond more closely to a 
realistic situation. 


Lateral Analysis 

The analysis of lateral control in the approach was per- 
formed for the status display and for the status- 
Interim-Director combination. In each instance allocation of 
attention was distributed according to the "optimal" values 
found in Chapters 3 and 4. That is, when the status display was 
being attended to, 75% of that attention was devoted to the 
localizer; when both displays were available, 80% of the atten- 
tion was devoted to the interim-director. 

Thresholds and residual noises were determined from [7]. 
Localizer thresholds were in anqular units, so that the effec- 
tive observation noise on lateral error and error-rate increased 
with distance from the runway- threshold . Thus, when approach 
was initiated, the lateral-error threshold was about 1.25m as 
compared to ,28m at the decision height. 

For the lateral approach analysis we assumed that the 
forward ground speed (i . e ., range-rate) was maintained at 31m./s 
This made the generation of the crosswind quite simple, i.e., 
we let 


v 

MW 


a h N = 


a U tan I\ 
o o 
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Because of the simplified dynamic representation used for the 
lateral analysis, the crosswind affected only the cross-track 
velocity of the aircraft. As opposed to the longitudinal case, 
no aerodynamic effects of the wind-shear considered. 

Pes u 1 ts . — "Window" performance for the two display con- 
ditions is shown in Table 24 and compared with correspond ina 
steady-state results. It should be noted that mean lateral errors 
were still decreasing but variances had very nearly converged when 
this "snap-shot" was taken. It can be seen that the mean-values for 
the shear cases are little affected by the display condition. 

Moreover, it turns out that the percentage change in variance 
from steady-state to wind-shear is also about the same for both 
display configurations. From an absolute standpoint, the director 
configuration shows considerably better performance in that smaller 
errors are achieved with less control. 

The above results are further illustrated in the time- 
histories of Ficrure 30. It can be seen that in the initial por- 
tion of the trajectory, when the wind is constant, the director 
configuration does better. After the shear starts the director 
continues to do better for some time and its peak mean lateral 
error is less than f or the status display. However, as 
the window is approached, the mean error becomes virtually ecrual 
for both displays. 

One might expect the director-system to show to advantaae 
with respect to workload renuirements . To test this notion, 
trajectories were obtained with a reference observation noise/signal 
ratio of -lldB for both the status display status and the status- 
director configuration. Scores at the decision-height for the two cases 
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FIGURE 30. Lateral Time Histories for STOLAND and Director 
Displays 
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Table 24 


LATERAL PERFORMANCE AT DECISION HEIGHT FOR 

CROSSWIND AND TURBULENCE 

(o = 1.5m/s, P = .20dB) 
v o 

a 



Status 

Status + 1 

Interim Director | 

Variable 

Steadv-State 

Crosswind 

EBIHSEB 

Crosswind 

y (m) 

0 

.73 | 

0 

.73 

0 (m) 

3.49 

4.13 

3.17 

3.77 

y 





y (m/s) 

0 

-.045 

0 

-.047 

(m/s) 

1.38 

1.61 

1.30 

1.47 

"J (deg) 

0 

-1.58 

0 

-1.59 

o , (dea) 

4) 

3.54 

4.36 

3.24 

3.90 

S’.Jdeg) 

w 

0 

-.40 

0 

-.4 

o. (dea) 

7.4 

9.2 

6.7 5 

8.27 

w 




• 
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(repeated for -20dB) are shown in Table 25. The STOLAND-Director 
combination is less sensitive to the change in noise/signal ratio. 
In terms of lateral error there is a smaller percentage increase 
in both the mean and standard deviation. Again, we see better 
performance for less control effort. 

Another aspect of director design that can be analyzed is 
beam capture. It w as decided to try a simple test of this aspect 
of the problem by introducing an initial (mean) displacement of 
10m from the localizer center-line. The small displacement was 
chosen, in keeping with the linearizing assumptions that were 
made. The results for status and Director cases are shown in 
Figure 31. The mean response for the two cases is indistingui- 
shable, so only one curve is drawn. However, the variability of 
the two cases is quite different as is seen by the plots of 
rms ( not standard deviation) values. For localizer intercept 
the director also appears to improve performance. 
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Table 25 

EFFECT OF NOISE RATIO ON APPROACH PERFORMANCE 

IN WIND-SHEAR 



Statv 

as 

Status-E 

>irector 

Variable 

P = -20dB 
o 

P = -lldB 
o 

P = -20dB 
o 


y 

.73 

1.73 

.73 

1.58. 

c y (m) 

4.13 

8.47 

3.77 

6.29 

y (m/s ) 

-.045 

-.09 

-.047 

-.08 

a* (m/s) 

y 

. 

1.61 

2.63 

1.47 

1.96 

$ (deg) 

-1.58 

-1.56 

-1.59 

-1.58 

a . (deg) 

<p 

4.36 

7.0 

3.90 

5.01 

6 w (deg) 

-.4 

-.39 

-.40 

-.39 

(deg) 

w 

9.2 

14.9 

8.21 

10.8 


SUMMARY 

Some aspects of the STOL approach in a mean-wind with shear- 
variation have been analyzed with the pilot-vehicle-display model. 
Results were obtained for longitudinal control with the status dis- 
play and for lateral control with the status display and the status 
interim-director combination. In general, the wind-shears degrade 
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performance by producing both mean errors and increased variabi- 
lity in the response, with the increased variability appearincr to 
be the major effect. 

Two wind profiles were considered in the analysis of longi- 
tudinal control, a decreasing and an increasina tailwind. Rela- 
tively good performance at the window was obtained in both cases. 
However, for the decreasina tailwind, the results showed that with 
the throttle fixed, excessive nozzle-control was renuired for wind 
compensation. When the nozzle control was limited (indirectlv, by 
penalizina control motions subseauent to shear-onset) , the height 
errors increased significantly. The rclativelv oood performance 
for the increasina tailwind \-ia c achieved wi th control reauirements 
that were not so excessive and, conseauontly , represent a more 
reliable result. The better performance is undoubted lv due to the 
additional time available for error compensation and is, of course, 
to be expected. 

The results for lateral control were about as expected. 

Specif icallv , the status-director combination performed better 
in terms of response to shears both from an absolute stand- 
point and in terms of sensitivity to observation noise/signal 
ratio. Beam capture response was also better for the configu- 
ration that included the director. 

To the extent that thev were invest iaated , the wind-shear 
responses tended to confirm essentially the steadv-state results 
(albeit that performance was worse in the wind-shear) . With 
regard to the details of the transient responses, it may be 
said that the behavior is easily explained. The transient res- 
ponse depends very much on the specific assumptions about initial 
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conditions and pilot strategy, which is not the least bit sur- 
prising. If one is interested in reproducing or predicting a 
particular time history (ensemble) then it is essential that 
conditions used in the model match those of the experiment. 
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CONCLUSION 

This chapter is comprised of two parts: for the convenience 

of the reader, we summarize the key results of the study and 
provide references to supporting data contained in the main body 
of the text; then we suggest areas for further work. 

Summary of Main Results 

A. With respect to the status displays of the STOLAND-EADI : 

(i) Steady-state analysis indicates that high levels of 
pilot workload (probably unacceptable) will be re- 
quired to achieve a 95% approach-success probability 
(Category II - window specifications) in median- 
turbulence conditions. When all turbulence levels 
are considered a 95% success probability is not 
possible with reasonable workload levels (Figure 13) . 

(ii) Even with stability auqmentation , the lateral-directiona 
control task is more difficult than longitudinal control 
To achieve a 95% probability of being within the lateral 
"window" requires approximately three times the "work- 
load" necessary to attain the same success- level for 
the height-airspeed window (Figures 6, 10). 

(iii) Comparison of proposed status displays with idealized 
status displays shows that the principal limitations 
imposed by the status display are the (a) the failure 
to provide adeouate information on alide-path and 
localizer error rates (i.e., vertical and horizontal 
flicrht path anqles) and (b) the requirement to share 
attention amona the display elements (Figure 8, Table 8) 
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(iv) Wind-shears, as expected, tend to degrade performance. 
They introduce mean-errors and increase response 
variability (Tables 22, 24). A decreasing tailwind 
with 2m/s / 30m (4KTS/100 1 ) shear- variation appears 
to exceed the rate-of-descent capability of the air- 
craft with throttle-fixed (Figures 27, 28). 

B. With respect to proposed display augmentation: 

(i) The proposed interim -flight directors provide 

improvement in performance at substantially reduced 
workload (Figures 17, 19, 20). When the average 
of all-winds is considered, the addition of the 
flight directors reduces the miss-probability by 
about a factor of two at all levels of pilot atten- 
tion (Figure 20) . To achieve a given probability of 
success for the median wind-condition, reauires 
about a third the workload with the directors as 
without them (Table 16) . 

(ii) Best performance on the longitudinal axis is obtained 
with attention shared roughly eoually between director 
and status displays (Table 11) . On the other hand, 
best lateral performance is obtained with an 80%-20% 
division of attention in favor of the director display 
(Table 15). Thus, the lateral director appears to 
provide a better "mix" of signals than does the longi- 
tudinal director, at least for the conditions analyzed. 
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( iii) The longitudinal director does not yield performance 
as good as that predicted by the idealized display 
(Table 12) , which fact suggests that the director 
laws could be further optimized. The lateral director, 
however, does lead to performance levels very close to 
those associated with the idealized display. 

(iv) Missed approaches with the directors, as without them, 
are primarily due to lateral errors (Fiaure 20) . In- 
sofar as the lateral director appears to yield per- 
formance close to that obtained with an idealized 
display, further improvement in performance via 
optimization of lateral displays will be very diffi- 
cult using the feedbacks currently proposed for the 
status and director displays. Although improvements 
in lonqitudinal performance will not increase approach 
success-probabilities appreciably, further reduction 
in longitudinal workload via display desian will allow 
more attention to the lateral task, thus increasing 
the overall success probabilitv or reducing the overall 
workload. 

(v) Response to cross-winds is improved by the lateral 

director, at least for the wind investigated (Table 24) 
The lateral director reduces the sensitivity of both 
mean and rms errors to attention-levels and will, 
therefore, reduce workload (Table 25) . 
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C. With regard to the pilot-vehicle-displav analysis procedure. 

(i) The "optimal-control model" for the pilot (with its 
interference and workload sub-models) in conjunction 
with modern systems analysis techniques, provides 
a powerful means for analyzing displays and their 
effect on system performance and reliability. 

(ii) The procedure may be used to determine bounds on 
expected display improvements via analysis of 
"idealized" displays (Table 5, 8). 

(iii) The model for task interference and workload permits 
linearized analysis of combined longitudinal and 
lateral performance in a rational and consistent 
manner. Interaction between axes is introduced via 
the limited capacity of the pilot and not through any 
vehicle coupling. 

(iv) The analysis techniques can serve as a basis for 

director design. Such a (preliminary) procedure for 
design of a longitudinal director system, that considered 
only the gust-regulation problem, yielded a configuration 
that resulted in substantially reduced workload 
(Figure 23 ) . 


Further Work 

As is customary, we conclude this report with suggestions 
for additional work and research. First, it is important that 
the results of this analysis be confirmed by simulation experi- 
ments. Beyond this validation, experimentation and data 
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analysis are needed to solidify the basis for choosing cost 
functional weightings when flinht directors are beina used. 
Additional approach data are needed to pin-down details of 
the pilot's time-varying adaptation. "Transient" data in the form 
of ensemble averages and variabilitv would be most helpful. 

Further display analysis seems warranted. For example, 
a more thorough model-analysis of wind-shear response with the 
director configurations is desirable. An interesting and 
relatively minor extension, that could not be accomplished within 
the constraints of this program, would be to take the estimation 
data from this study and analyze pilot decision-making with 
regard to missed approaches [12] . Display performance in curved 
approaches is also an important area for investigation that 
is amenable to the technioues employed here. 

It appears that the approach to flight-director design 
suggested in Section 5 has much merit and is worthy of continued 
investigation. It became apparent in that analysis that the 
facility for ignoring certain model feedbacks would be useful 
and the programs should be modified accordingly. 

Finally, it is important to take cognizance of the changing 
role of the pilot from one of active controller to monitor and 
manager. Analysis procedures such as those employed here are 
necessary for investigating and understanding this role. It is 
our conviction that the models and techniaues used here can serve 
as a basic building block for the development of models for flight 
management and that such development should be addressed vigorously. 
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APPENDIX A 

Mod ification for Tiire-Varving Input Disturbances 

In Reference 6, procedures for modifying our model of the 
human operator to account for constant mean disturbance inputs 
were developed. Here, we consider a wider class of inputs, viz 
those that have a time-varying mean. Specifically, we consider 
the class of disturbances that may be modelled as the state (or 
output) of a linear dynamic system: 

z = ; £(t Q ) = z c ( A . 1 ) 

where A , a constant, and z are assumed known*. Various distur- 
— z - o 

bance inputs may be generated from such a model. For a detailed 
discussion, see Ref. 19. The results we present appear in various 
forms in the literature [Ref. 19, 20] and are developed here for 
completeness and convenience. 

The disturbance state is assumed to enter the system 
linearly, so the system state eouations are 

x=/vx+Bu + Ew + Fz_ (A. 2) 

The problem is to choose u to minimize 

J = r;{i /[x* 0 X + u' R u]dt} (A. 3) 

1 1 o 

where Q > 0, R > 0. 


+ Thus , the state z (t) is known for all t. In applying the model, 
jz(t) v/ill have to be estimated from available, noisy outputs, as 
are other system states. 


A-l 
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To solve this problem, we define an augmented state 

x' = [z'ix’l satisfying 
I 


x = 


z_ 

X 


A O z 
—z J_— _ _ r _ 

F I A ! ; x 


0 j 0 

+ u + w 


■ B 


e ; 


(A. 4a) 


A x + B u + w , 


(A. 4b) 


an auamented weicrhtina matrix 


r o i o 1 



0 I o 


and consider the minimization of (A. 3) for fixed T. This is a 
well-def ined problem whose solution is given by 


U*(t) = -R 1 £' K ( t ) X(t) 

where 

-K = A' K+KA + Q- KB R _1 B' K ; t<^T 
K (T) = 0 

Letting K be partitioned conformally with A, viz. 


■hi ! 

-_-21 

L — 21 1 

—2 2 . 


(A. 5) 


(A. 6) 


(A. 7) 


A- 2 
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we find that 


u*(t) 



(K 21 (t)z + K 22 (t)x) 


(A. 8) 


and 



-11 


A 


A' 

— z 


*11 


+ —21 


F + F 


- 21 ' 


*21 


B R 


-1 


B ' K 


21 


K U (T) = 0 


(A . 9a) 


■hi ■ hi ±z + h2 E + - hi ~ -? -1 E ' )K 21 ; 


k 21 ( T) = 0 


(A. 9b) 


'*2 2 


*22 t + K 22+ Q 


*22 E E' 1 E' E„ . 


K. 22 (T) = o (A. 9c) 

Mote that Fan. (A. 9c) is the P.iccati Fan. for the regulator 
problem without the time-varvine mean disturbance. Also, since 
(A. 9b, c) and (A. 8) are independent of the solution to (A. 9a) 

is not needed in implementing the optimal control. 


The structure of the controller implied by (A. 8) and (A. 9) 
is shown in Figure Al. The feedback structure of the controller 
is identical to that of the optimal linear regulator; the modifi- 
cation to the regulator is a set of feedforward gains operating 
on the disturbance (on the estimate of the disturbance, in practice) 

We wish to explore conditions that lead to constant gains 
in (A. 8) inasmuch as this simplifies the model implementation and 
interpretation considerablv . Thus, we consider the case, T -*•<». 


A-3 
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It is well-known that the solution to (A. 9c) yields 


lim K_,,(t) = K = constant (A. 10a) 

T-*-oo 22 


where 


K A + A' K + Q 


K B R 


E 'K 


0 


(A. 10b) 


Using (A. 10), it can be shown that 

lim K„, (t) = K = constant (A. 11a) 

T-m» 2 

satisfying 

K A + A' K = -K F 
— z — z — — z — — 

A = A - B R -1 B' K 

if, and only if. 

Re {A^ + y_.) <0 V i, j (A. 13) 

where A^ and y^ are the eigenvalues of A and , respectively. 
Inasmuch as Ee (hi <0 V i by virtue of the stabilizing properties 
of the optimal regulator, a sufficient condition for (A. 13) to 
hold is that 

Re (y.) < O , V j (A. 14) 

The condition that the disturbance not have an exponentially 
growing component (which is the case if Re < O) is sufficient 

for the gains on the disturbance to be constant, but does not 


(A. lib) 
(A. 12) 
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guarantee the optimality of our solution. This is so because 
the cost functional may not be finite as T-*“. This would not 
seem to be a serious complication for our purposes. The upper 
limit T may be considered to be sufficiently large for the solu- 
tion to the Riccati equation to converge: Once the gains have 

been computed/ as constants , we evaluate the various scores by 
integrating appropriate equations forward in time. For realistic 
inputs and reasonable times, the scores will be reasonable. Of 
course, if the cost functional (A. 3) is bounded as T-*°°, then we 
have the optimal solution. 

In the special case of a. constant disturbance, A =0, and 
Eqn. (A. 14b) yields 

K z = "(A)" 1 K F (A. 15a) 

and, substituting for u in (A. 2) and setting W E 0, we find 

lim. x ( t ) = - A^F - B R” 1 B’ K ) z = constant (A. 15b) 

ts-oo ~ “ “ ~ ~ z ~° 

Equations (A. 15a) and (A. 15b) verify the results obtained "directly" 
in Ref. 1. 
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APPENDIX B 


Computation of System Response to Specif i c Initial Condit ions 

In this section we develop general expressions for computing 
model responses to specific initial conditions"*". These formulae 
are needed to compute control and monitoring performance associated 
with specific time-varying disturbances, e.g., mean winds (see 
Chapter 6). For simplicitv, we consider here the case where 
the human's time-delay is negligible (or may be accounted for by 
adjusting t. t ) , i.e., x = 0. 

We begin by assuming that, in general, the initial state 

x(t ) is a gaussian random variable with mean m and covariance 
-o'- - o 

X Q . The equations governing the motion of the optimal closed 
loop system: are well known and given by [6]. 

A 

x = £ 2 - B L* x + E w(t) (B. la) 

- 

x = F x + P (t) [C (t) e + v ( t ) ] ; x(t ) = m o (B. lb) 

e = C(t) e + E w(t) - P(t) v (t) (B.lc) 

v’here 

F £ A - B L* 

G(t) A A - E(t)C(t) V“ 1 (t) C(t) ( B . 2 ) 

P (t) A Z (t) C’ (t) V _1 (t) 


+ These results are a generalization and amplification of those 
given in Appendix C of Pef . 6. In addition, the derivation is 
somewhat different and, hopefully, more straightforward. 


B-l 
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and 


+ 

f 


W(t) 6(t-x) = (w(t)w T (x)} 


V(t) 6 ( t-x ) = {v(t) v T (t)} 


(B. 3) 


L* is the "optimal gain" matrix computed from the optimal 
linear regulator problem and £(t) is the solution to the variance 
equation. 


£ = A £ + £ A'- £ C' (t) V (t) C ( t) £ + E W(t) E' ; 


1 (t ) = X 
— o — o 


which may also be written as 


> (B . 4 ) 


£ = G ( t) £ + £ ( t) G' + P(t) V(t) P’(t) + E V7 ( t) E' 


= *o 




In most problems it is assumed that V(t) and W(t) are known 
a priori, so £_(t) may be precomputed. In the case of the human- 
operator-model application, V(t) depends on the state and £Jt) 
is computed on-line; this is possible because solution to (B.4) 
is carried out in the forward direction. In either instance, the 
optimal filter implementation involves solution to (B.l) and (B.4). 


For convenience, we define the composite "state" and "input as 


+ 


Bar denotes expectation operator. 


B-2 
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X = 

X 

/N 

X 

; w = 

- ~ 

w 




V 

Nv 


Then, (B. 1) may be written as 
x = A { t ) x + E ( t) w 

where 


A 

-BL* 

0 


E 0 

0 

F 

P(t) C(t) 

; E = 

0 P(t) 

0 

0 

G(t) 


E -P(t) 

L 






We now consider the response of (B.6) to a s pecif 
condition x(t Q ) drawn from the original distribution. 


£ (t o> 


m. + e 
-o — o 


Thus , 


f m + e 

— o — o 


m 

— o 


L 


Then 


_ t 

X ( t ) = S' (t , t ) X ( t ) + f ,y (t,T) E (T) W ( T ) dT 

- - 0-0 t ~ “ • - 

O 


- x(t o } = 


( B. 5) 


(B.6) 


( B. 7 ) 


ic initial 


( B. 8 ) 


( B. 9 ) 


(B.10) 
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where 


d - (t,t Q ) = A(t) l (t,t o ) ; ¥(t 0 ,t 0 ) = I 

dt 

Writing H* in partitioned form corresponding to A, 


y 

-11 

y 

-12 

y 

-13 

y 

-21 

y 

-22 

y 

-2 3 

y 

-31 

y 

-32 

y 

-3 3 


(B. 11) 


it is relatively straightforward to compute the sub-matrices 

V . . , in terms of y__ owing to the trianqular form, of A. The 
-i] —33 — 

results are presented in Table Bl. 


We are actually interested in the ensemble of responses, 
to this specific initial condition, that are generated by 
samples of the random processes w( . ) and v(.). These are 
completely characterized, because of our assumptions, by the 
first- and second-order moments, i.e., the mean and variance. 


The mean-response of (B.6) given (B.9) is readily calculated 
from; (B.10) 


x(t) = ¥(t,t Q ) x ( t o ) + | V(t,T) E(t) w(T)dx 


If w(.) = 0, as we will assume, 


x(t) = V(t,t Q ) x(t Q ) 


(B. 12) 


We note that starting from; a specific initial condition x(t ) = 
x(t Q ) . We have included the expectation on initial conditions for 
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TABLE B1 


d^ 

— 33 ^t' x ) : Solution of — (t,x) = G(t) '*' 33 (^ 1 )#’ 


dt 


-33 


—3 2 ( t » T ) = 0 


2 31 (t,T) = o 


2 23 (t,T) = / e- (t S) P(s) C (s) j ^ 33 (s,i)ds 


T 22 (t,T) = e 


F(t-x) 


I 21 (t,t) = 0 




S' (t , 1 ) = - f e- (t S) BL* e- (s T 5 ds 
12 T 


t s 

V (t,T) = - / ds e — (t-s) BL* [}e- (s_0) P(c) C(o) ^ 33 (c 

13 T T 


,T) = I 


,T ) do ] 
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generality. In cases where there may be step changes in distur- 
bances, t could correspond to a time after the initial time and 
x(t ) would be a random variable. A siir.ilar consideration will 
be necessary in computing the variances. 


In order to compute the covariance of x(.), we define 


X = 


— — 


— 

£ 


x - X 

n 

= x - X = 

a 7T 

x - X 

Y 


e - e 



mm~ — 


( B. 13 ) 


It is easily seen that 

X = A ( t) X + E ( t ) w 


T — A 

and, moreover, that E { X_ X_ } = cov {x} = A, satisfies 


dA = 
dt 


A ( t ) A + A A' (t) + E ( t ) fi(t) E'(t) ; A (t )= A (fi.14) 

— — — — — — — — o —o 


where 


fl(t) = 


W(t) | 0 

0 1 V(tJ 


( B. 15 ) 


The solution to (B. 14) may be written directly as 

A(t) = I(t,t o ) A(t 0 ) ¥' (t,t Q ) +Jv_( t.r) E(x) £(t) E ' (t ) 

fc o 

V 1 (t ,r ) dr 


( B. 16 ) 


In practice, because of the triangular form of A, it is not 
necessary to obtain all sub-matrices of A_ (t) to determine cov {x} , 

/v 

cov {x}, and cov {e}; indeed, the only additional computation 
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involves determining cov {x, e } and its transpose. Thus, 
(B.14) yields 

cL cov {e} = d _-33 = G(t) A 33 + A 33 G 1 (t) + E W(t) E' + 
dt dt 

P(t) V(t) P' (t) 

d cov { x , e }= d — 2 3 = F A_ 2 3 + A 23 — ' (t) + P(t) C(t) 
dt “ dt 

-P(t) V(t) P'(t) 

A 

A = cov {e, x) = A_2 3 

d cov {x} = d --2 2 = F A 3 2 + A 22 A' + P(t) C(t) A32 
dt dt 

+ [P(t) C(t) A 32 (t)]' + P(t) V(t) P' (t) 

d cov {x} = ^lj. = A A 33 + A A’“[B L* ^22^^ 

dt d t ‘ 1 

- [B L* A 22 (t) )' + E W(t) E’ 


(B. 17) 


(B. 18a) 


(B. 18b) 


( B. 19 ) 


(B.20) 


Comparing Err ns . (B. 17)and (B.4), we find 

d_ (A 33 - E) = C ( t ) (A 33 - E) + (A 33 - E) G' (t) 

dt 

thus 

A 33 (t) = cov {e (t) } = E ( t ) + i(t,t Q ) [A 33 (to) - Zit Q )] 

$' (t ,t Q ) (B. 21) 
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where 


!<t,t 0 ) = G(t) £ (t,t D ) ; £ (t o ,t o ) = I 


(B. 22 ) 


Eqns. (B.18) may be integrated directly to obtain 

t 


A 2 3 ( t ) = e- (t V A 23 (t D ) £’ (t ,t Q ) + 


e^ (t " T) P(T) 


ECU) A 33 (t) - V ( T ) P' ( T ) ] £' ( t ,T ) dt 


or , using ( B. 2 ) , 


A 23 (t) = A 23 (t Q ) £ , (t,t o ) + f e— (t x) - (t) £ (t) 


[£ 33 (t) “ 1 1 ' (t ,T)dx 


(B.23) 


Substitution of ( B. 2 1) into (B.23), gives 

cov{x(t) ,e(t) } = A 2 3 (t) = e ~ (t t °' > A 2 3 + 


^ 32 (t Q ) ~ A(t o )l £Mt,t o ) 


( B. 2 4 ) 


where 


£ ( t , t Q ) = / e — ^ P(t) C (t ) £(T,t Q )dT 


(B. 25) 


. + 


In a similar fashion, we obtain 


Alternative expressions are, of course, possible. For example 
one can^solve for cov {x} = cov (x) + cov {e} + cov {x,e} + 
cov {e,x}, as was done in Ref. 6 . 


B -8 
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cov {x(t) } = A 22 (t) = e— ^ t o ) ^ 22 ^ 0 ^ e ~ ^ to ^ + 


ana 


J e- (t_T) [P ( T ) C ( t ) A 32 (t) + A_ 2 3 ( x ) C'(t).(B. 26 ) 
t 0 ( 

P 1 ( T ) + P ( T ) V ( T ) P ' ( T ) ] e F ' (t_T) dT 


cov {x (t) } 


= e- (t t o ) A, , (t ) e-- (t t o ) + 
—11 o 


e A(t-T) 


[E W(t) E' - (B L* A 22 (t) ) 


( B. 2 7 ) 


A * ( f" — T ) 

(B L*A 22 (t) ) *] e- v T; dt 


Eqns. (B.21) and (B.24) are especially revealing. We see 
from (B. 21) that if 


A,, (t ) = I (t ) = X ( B. 28) 

— 3 3 o o —o 

then LJt) is the actual covariance. Also, referring to (B.24), 

if A_ 2 3 ( t Q ) = 0 and (B. 28) holds then the estimate and the error 

are uncorrelated for all t. When x(t ) is considered as a random 

— o 

variable the above conditions hold and the filter is optimal with 
respect to the given data. 

For a sample path corresponding to a specific initial 
condition (or for the ensemble of paths corresponding to that 
condition), ^ 33 ^ 0 ^ = 0 hecause e(t- 0 ) f; “- xe< 3- Thus, (B.28) does 

not generally hold, the actual error covariance differs from £ and 
the error and the estimate are correlated. The filter is optimal 
in the same sense as before, i.e., based on the prior knowledge 
of the random variable x(t Q ). However, one could of course do 
better, if the initial condition were known exactly. 
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The above results are interestincj from another standpoint. 
It is possible that A 33 (t 0 ) ¥ £(t Q ) even though we are not 
considering a specific initial condition. Such would be the 
case, e.g., if the pilot did not know the initial state covariance 
perfectly, as is most certainly the case. Enuation (B.21) then 
reveals the sensitivity of the solution to this type of error. 
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APPENDIX C 
Vehicle Dynamics 

In this Appendix, we present longitudinal and lateral dyna- 
mics that were used in the analyses described in the report. The 
basic equations were linearized perturbation equations of standard 
form [21] . 


Longitudinal Dynamics 


In vector-matrix notation, the longitudinal dynamics are 
x = Ax + B_6 + Ew 
y = C x + d 6_ 

/ 

The state, control and output vectors for longitudinal control 

★ 

were for the no-director case 



- = (u g' w g' u ' w ' G ' °' h) 


A = <6 e , V 


y = (h, h, e, ct, u) 
where superscript T denotes transpose, 


/ 


C-2 


Two sets of longitudinal dynamics were investigated in this 
report because of a change in stability derivatives and trim- 
conditions that was made in the course of the effort. Initial 
values for the system matrices as provided by Ames personnel 

*When a director is used that has dynamics (from filterinci, 
e.g.) additional states are needed. The additions are not 
shown here. 
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were used for the analysis of the STOLAND DISPLAY (ALONE) and 
the '"model-based" director. These dynamics are given by the 
computer printout labeled STL3 (pp4-5) .The second set of dyna- 
mics were used to analyze the "interim" director (because the 
interim director was designed for these dynamics. These dyna- 
mics are given by the computer printout labeled STL9. 

It is of interest to compare performance for the two sets 
of dynamics. This was done for the 1%-wind using nominal condi- 
tions (the STOLAND display; the resulting rms performance scores 
pgy*$a in fahip pi , Thp f f pv-pppps are clearl^ ^inor. 

Table Cl 

COMPARISON OF PERFORMANCE SCORES FOR TWO 
SETS OF LONGITUDINAL DYNAMICS 



C-2 
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Lateral Dynamics 


Simplified lateral dynamics illustrated in the block 
diagram of Figure 18, were used. The lateral states; controls, 
and outputs are 

X T = (V q , p, ({), y) 


• • * 

y = (y» y> 4>r 4>) 

The values for A, B, C, E in the lateral case are given by 
LAT1 for the SAS-off case (p. 8) and LAT3 for the SAS-on dyna- 
mics (p. 9 ) . 


*The vector y is used to denote, generically, 
whereas the scalar y is the lateral error. 


the display vector 



Report No. 2484 


Bolt Beranek and Newman Inc 


FILE NAME: STL3 


TOTAL NO. OF 

NOISE STATES= 

2 



A MATRIX: 





- 1 . 500E-0 1 

0.000E-01 

0.000E-01 

0. 000 E~ 0 1 

0.000E-01 

0.000E-01 

0.000E-01 




0 » 000E-0 1 

“ 1 * 000E+ 00 

0. 000E-01 

0.000E-01 

0.000E-01 

0.000E-01 

0.000E-01 




-3* 600E-02 

1 . 600E-01 

-3. 600E-02 

1 . 600E-01 

“ 1 .700E-01 

0 * 000E- 0 1 

0.000E-01 




“2* 870 E- 01 

-4.860E-01 

-2.870E-01 

-4.860E-01 

1 .030E-02 

5.320E-01 

0.000E-01 




0.000E-01 

0. 000E-0 1 

0.000E-01 

0.000E-01 

0.000E-01 

1 . 000E+00 

0.000E-01 




4.360E-01 

-2* 230E-0 1 

4. 360E-01 

“2*230 E- 01 

-7. 580E-03 

- 1 .2R0E+00 

0 » 000E- 0 1 




0» 000E- 0 1 

0 • 000E- 0 1 

-6. 120E-02 

“ 1 » 000E+ 00 

5- 340E-01 

0.000E-01 

0 • 000E- 0 1 





B MATRIX: 

0. 000E-01 
0. 000E-01 
0. 000E-01 
-2. 390E-02 
0 • 000 E- 0 1 
- 1 . 120E+00 
0. 000E-01 


E MATRIX: 

5» 500E-0 1 
0 • 000 E- 0 1 
0.000E-01 
0.000E-01 
0 • 000 E- 0 1 
0. 000E- 0 1 
0.000E-01 


0 • 000E- 0 1 
0. 000E-01 
-2. 860E-02 
6. 670E-04 
0* 000E-0 1 
-9. 1 70E-02 
0. 000E-01 


0.000E-01 
1 . 4 J 4E+00 
0» 000E- 0 1 
0. 000E-0 1 
0. 000E-01 
0. 000E-01 
0. 000E-01 


04 



Report No. 

2484 
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0. 000E-01 

0. 000E-01 

0. 000E-01 

0 • 000E- 0 1 

1 • 000E+ 00 


0. 000E-0 1 

0. 000E-0 1 

-6. 100E-02 

0. 000E- 01 

0. 000E-01 


0 • 000 E- 0 1 

0 • 000E- 0 1 

0. 000E-01 

0 • 000E- 0 1 

0 • 000E- 0 j 


0.000E-01 

0. 000E-0 1 

0. 000E-01 

1 • 000E+ 00 

0. 000E-01 


1 • 000 E+ 0 0 

0 • 000E~ 0 1 

1 • 000E+ 0 0 

0 • 000E- 0 1 

0. 000E-01 



Bolt Beranek and Newman Inc. 

0.000E-01 0.000E-01 
1 • 000E+ 00 5.330E-01 
0.000E-01 1.000E+00 
0.000E-01 0» 000 E- 01 
0«000E-01 0»000E~01 


D MATRIX: 

0. 000E-0 1 
0. 000E-0 1 
0.000E-01 
0.000E-01 
0. 000E- 0 ] 


0. 000E-01 
0. 000E-01 
0. 000E- 0 1 
0. 000E-0 I 
0. 000E-01 
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FILE NAME: STL9 


TOTAL NO. OF NOISE STATES 3 2 


A MATRIX: 


- 1 . 500E-01 
0.000E-01 

0. 000E-01 
0.000E-01 

0. 000E-0 1 

0.000E-01 

0. 000E-01 

0.000E-01 
0 . 000 E- 0 1 

~ 1 o 000E+ 00 
0.000E-01 

0.000E-01 

0.000E-01 

0 * 000E~ 0 1 

-5.200E-02 
0 * 000E" 0 1 

1 .230E-01 
0.000E-01 

-5.200E-02 

1 .230E-01 

- 1 . 700E-01 

-2. 760E-01 
5. 300E-01 

- 5.220E-01 
0 * 000 E~ 0 1 

-2.760E-01 

-5.220E-01 

1 . 740E-02 

0. 000E-01 
1 . 000E+ 00 

0.000E-01 
0. 000E-01 

0.000E-01 

0.000E-01 

0.000 E- 01 

5.020E-01 
- 1 . 350E+00 

-2* 040E-0 1 
0.000E-01 

5.020E-01 

-2 • 040E- 0 1 

- 1 • 200E- 02 

0.000E-01 
0 * 000 E~ 0 1 

0.000E-01 
0. 000E-01 

2.270E-02 

" 1 » 000E+ 00 

5-370E-01 


B MATRIX: 

0 . 000 E~ 0 1 
0 . 0(^0 E“ 0 1 
0. 000E-01 
-2. 400E-02 
0. 000E~ 0 1 
- 1 . 280E+ 00 
0.000E-01 


E MATRIX: 

5. 500E-01 
0. 000E-01 
0. 000E-01 
0 • 000 E- 0 1 
0» 000E** 0 1 
0 * 000 E“ 0 1 
0. 000 E- 0 1 


0. 000E-01 
0. 000E-01 
-2.8R0E-02 
2»220 E- 03 
0.000E-01 
- 5.630E-02 
0 * 000E- 0 1 


0 * 000E~ 0 1 
1 . 41 4E+00 
0 . 000 £“ 0 1 
0 » 000E - 0 1 
0 * 000E- 0 1 
0.000E-01 
0 » 000E- 0 1 
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C MATRIX: 


0. E- Pt 1 

0.000E-01 

0. 000E-01 
1 • 000E+ 00 

0. 000E-01 

0.000E-01 

0.000E-01 

0. 000E-0 1 
0.000E-01 

0.000E-01 
0. 000E-01 

2* 270E-02 

- 1 . 000E+00 

5.370E-01 

0.000E-01 

0.000E-01 

0.000E-01 
0. 000E-01 

0. 000E-01 

0. 000E-01 

1 • 000E+ 00 

0.000E-0 1 
1*000 E+ 0 0 

0. 000E-01 
0. 000E-01 

0. 000E-01 

0. 000E-01 

0. 000E-01 

1 . 000E+00 
0.000E-01 

0.000E-01 
0. 000E-01 

1 • 000E + 00 

0.000E-0 1 

0.000E-01 

MATRIX: 

0 • 000 E- 0 1 
0* 000E- 0 1 
0.000E-0 1 
0.000E-01 
0. 000E-0 1 

0. 000E-01 
0. 000E-01 
0. 000E-01 
0. 000E-01 
0. 000E-01 
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Bolt Beranek and Newman Inc 


FILE NAME: LAT 1 


TOTAL NO. OF NOISE STATES= 1 
A MATRIX: 

-1.500E-01 0. 300E-0 1 0.000E-01 0.300E-01 0.000E-01 

8.540E-02 -5.964E-01 2.500E-01 0.030E-01 0.000E-01 

0.000E-01 1.000E+00 0.000E-01 0.000E-01 0.000E-31 

0.000E-01 0.000E-01 3.160E-01 0.000E-01 0.000E-01 

0.000E-01 3.000E-01 0.0 00 £"01 5.410E-01 0.000E-01 


B MATRIX: 

0. 000E-01 
X. 1 1 5E-01 
0.000E-01 
0. 000E-01 
0.000E-01 


E MATRIX: 

5. 500E-01 

0. 000E-01 
0.003E-01 
0.000E-01 
0 • 0 3 0 E - 0 1 

C MATRIX: 

3.003E-01 0.000E-01 0.000E-01 0.000E-01 I • 000E+00 

3 • 000E-0 1 0.000E-01 0.000E-01 5 • 4 1 0 E - 0 1 0.000E-01 

0.003E-01 0. 003E-01 1.000E+00 0.000E-01 0.000E-01 

0.000E-01 1 • 300E+00 0.008E-01 0.030E-01 3.000E-01 


0 MATRIX: 

0. 030E-01 
0.000E-01 
0* 000 E “01 
0.000E-01 

) 
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Bolt Beranek and Newman Inc 


FILS NAME: LAT3 


TOTAL NO. OF NOISE S T A T E S = 1 


A MATRIX: 

- 1 • 530F.-0 1 
0. 000 E -01 
0 . 003E-0 1 
A . A 91 A E - 0 1 
1 . 303E+33 


8 MATRIX: 

0.000E-01 
6. ROME -01 
0.000E-01 
0. 300E-01 
0. 000 E- 01 

E MATRIX: 

5. 509E-01 
0 . 000F.-0 1 
O • 000E-0 1 
0.033E-01 
0 • 000 E _ 0 1 

C MATRIX: 

0. 000E-01 
1 • 000 E +00 
0. 000E-0 1 
0 • 0 0 0 E — 0 1 


0.000 E-S 1 
-1 . 68 3 E +00 
1 . 300E+0W 
3.003F.-31 
0 . 030 £-13 1 


0. 000E-01 
0. 330E-01 
0.000E-0 1 
1 .00 3 E + 3 id 


0.003E-01 
0 .000E-0 I 
-9. 600 E -08 
3.210E-01 
0.030E-01 


0 • 0 0 3 E - !'} 1 
0.033E-3 1 
1 . 033 £+33 
-9. 603 E- 02 


3 * 0 0 0 E - 3 1 
0 . 0 0 0 E - 0 1 
0 • 303E-0 1 
0.303E-Q 1 
5. 3 ROE- 01 


0.000 £-31 
5. 323E-0 1 
0.000 £-01 
3 • 0 3 0 E - 3 1 


s ) MA T k I X : 

0. 3 3 3 E - 0 1 
3 • 3 0 0 F - 3 1 
0.003E-01 
0. 000E-31 


0 * o y o e - 0 1 
•3 • 303 £-0 1 
3. 000 £-3 1 
0 . 0306.-0 1 
0.000E-0 1 


1 . 003 £+00 
0 . 0 0 O E - 0 1 
0 . 000 £-01 
0 * 000 E“0 1 
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